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ABSTRACT 

The  radiational  model  employed  here  computes  the  planetary 
albedo  and  the  solar  absorption  by  atmospheric  layers  and  by 
earth's  surface  for  the  primitive  equation  model  of  FNWC. 
Large-scale  cloud  parameterization  in  several  layers  enters 
into  these  computations.   The  solar  computations  are  made  on 
a  gridpoint  basis  employing  the  water-vapor  mass  over  each 
gridpoint.   Longwave  cooling  effects  over  the  same  layers  and 
at  the  earth's  surface  are  calculated  using  emissivity  formulas 
after  Sasamori,  and  are  also  dependent  upon  the  cloud- 
parameterizations . 

Two  forms  of  the  cloud-parameterizations  were  tested  using 
FNWC  gridpoint  oceanic  data  for  16  January  1974.   The  objective 
was  to  determine  the  parameterization  which  better  verified 
the  radiational  balance  as  a  function  of  latitude  at  the 
tropopause,  as  compared  with  the  January  satellite  climatology. 
Separate  heat-budget  checks  of  the  tropospheric  columns  and  of 
the  oceanic-water  mass  were  tested  for  consistency  with  January 
climatology.   The  smaller  cloud-parameterization  values  gave 
the  best  verification. 
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I.   INTRODUCTION 

This  thesis  represents  a  continuing  study  of  a  radiative 
heating  parameterization  for  use  in  appropriate  a-layers  of 
the  Fleet  Numerical  Weather  Central  (FNWC)  numerical  predic- 
tion system.   The  immediate  goal  of  this  thesis  is  to  analyze 
the  heat  balance  of  the  earth-ocean  system  utilizing  FNWC 
gridded  data-fields  at  constant  pressure  levels  for  16 
January  1974.   These  gridpoint  data  were  selected  along  four 
meridians,  three  in  the  Pacific  and  one  in  the  Atlantic 
Ocean.   Most  of  the  gridpoints  were  located  in  the  Northern 
Hemisphere  (Fig.  1) . 

The  most  important  influence  on  the  radiational  disposi- 
tions (short  and  long-wave)  of  this  study  was  the  specifica- 
tion of  amounts  of  clouds  in  two  specific  layers.   The 
initial  specification  of  the  fractional  amounts  CL(1)  and 
CL(2)  are  based  on  large-scale  formulations  developed  by 
Smagorinsky  [1960]  and  first  used  in  a  similar  study  based 
upon  16  October  1973  data  by  Warner  [1974]. 

The  governing  equations  of  the  radiational  transfers  of 
the  model  have  been  devised  by  Martin  [1972,  1974],  who 
modified  the  streams  of  radiation,  both  solar  and  terrestrial, 
to  respond  to  the  presence  of  clouds  in  specified  amounts  in 
the  two  layers.   The  radiational  model  developed  by  Martin 
and  utilized  here  has  similarities  to  those  in  use  in  UCLA 
and  NCAR  General  Circulation  models,  but  with  special  adapta- 
tion to  the  a-levels  of  interest  to  FNWC. 
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Figure  1.   FNV/C  polar  stereographic  grid  and  meridians  (lines 
1,  2,  3,  and  4)  selected  for  study.   The  longitudes  A  are 
shown  for  each  meridian  as  well  as  the  extent  considered  of 
each  meridian. 
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Warner  found  that  his  global  albedo  estimates  were  too 
high  as  compared  with  the  satellite  climatology  of  Raschke 
et   al .     [1973],  thus  leading  to  underestimates  of  the  radia- 
tive balance  both  at  the  top  of  the  atmosphere  and  at  the 
surface  for  the  data  of  his  experiment.   In  this  study,  the 
original  Smagorinsky  [1960]  cloud-specifications  CL i  and  CL2 
were  computed  for  use  in  the  model;  however,  reduced  cloud- 
amounts  (2/3-CL)  in  each  layer  were  also  included  for  radia- 
tive budget  comparisons. 

The  comparative  results  afforded  by  the  "2/3-CL"  model 
gave  reasonably  close  agreement  as  a  function  of  latitude 
with  the  radiative  climatology  of  Raschke  et   al.     [1973] , 
for  the  NIMBUS  III  period  21  January  -  3  February  1970, 
examined  by  these  experimenters. 

In  order  to  analyze  the  effects  of  the  "2/3-CL"  versus 
the  "full-CL"  model  used  by  Warner,  oceanic  and  tropospheric- 
heat  balances  were  also  computed  for  both  sets  of  cloud  con- 
ditions for  16  January  1974.   It  was  necessary  in  this 
connection  to  augment  the  purely  radiative  model  by  a  turbu- 
lent boundary-layer  model  for  sensible  and  latent  heat 
transports  as  adapted  from  Kaitala  [1974] . 

In  the  2/3-CL  case,  the  zonally  averaged  results  for 
the  components  of  the  ocean- tropospheric  system  gave  realis- 
tic results  when  compared  to  climatology  for  this  mid-winter 
period. 

The  possible  error-source  of  using  one  set  of  computation 
al  data  to  be  representative  of  mid-January  climatology  is 
recognized.   Nevertheless,  the  heat-package  results  computed 
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with  the  2/3-CL  cloud  model  were  encouraging  when  examined 
against  climatology.   The  major  effect  noted  by  the  reduction 
in  cloud  cover  was  the  reduction  in  global  albedo  in  both 
low-  and  mid-latitudes  thus  contributing  to  greater  surface 
net  heating  rates.   This  latter  result  is  in  general  accord 
with  recent  observations  of  Von  der  Haar  and  Hanson  [1969], 
Von  der  Haar  and  Oort  [1973],  as  well  as  with  Raschke  et    at. 
[1973].   However,  the  specification  of  layered  amounts  of 
clouds  in  the  tropics  still  seems  to  give  somewhat  excessive 
global  albedos.   It  is  felt  that  the  geometry  of  solar  radia- 
tion streams  impinging  upon  cumulus  subgrid  cells  in  the 
tropics  must  be  reexamined  so  as  to  divert  a  greater  per- 
centage of  solar  radiation  downward,  as  compared  with  cor- 
responding cloud- layering  reflective  effects  in  middle  and 
high  latitudes. 
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II.   DATA  PREPARATION 

Temperature  and  humidity  data  utilized  in  this  study 
were  taken  from  gridpoint  soundings  along  four  meridians 
of  the  FNWC  63-by-63  gridmesh.   The  data  gridpoints  along 
these  meridians  (Fig.  1)  were  chosen  so  as  to  ensure  condi- 
tions representative  of  oceanic  locations  for  computations 
of  the  radiational  and  other  heat  budget  items.   The  reasons 
for  choosing  oceanic  locations  for  such  computations  were 
threefold:   (1)  the  maritime-area  heating-rate  computations 
are  likely  to  be  relatively  stable  timewise  over  the  data- 
day  (16  January  1974);  (2)  the  constant  a  =  P/tt  surfaces 
of  the  FNWC  primitive  equation  system  are  close  to  being 
constant  pressure  levels;  and  (3)  values  of  the  heat-budget 
items  computed  here  may  be  compared  with  corresponding  cli- 
matology during  mid-winter.   Corresponding  values  over  land 
areas  appear  to  be  less  amenable  to  comparison. 

Along  the  four  oceanic  meridians  used,  a  total  of  67 
gridpoint  soundings  were  selected  from  I  =  1  to  I  =  25  on 
lines  1,  2,  and  3  (over  the  Pacific  Ocean),  and  26  soundings 
from  I  =  63  to  I  =  38  along  the  Atlantic  meridian,  line  4 
(cf . ,  Fig.  1).   The  device  of  selecting  gridpoint  soundings 
chosen  from  the  indicated  meridians  of  the  FNWC  polar 
stereographic  map,  made  it  unnecessary  to  resort  to  spatial 
interpolation  between  original  data  gridpoints  along  the 
meridians . 

The  gridpoint  soundings  from  the  data  lines  just  mentioned 
were  taken  from  the  original  FNWC  63-by-63  gridpoint  analyses 
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of  temperature  at  the  surface  and  of  T(P)  and  q(p)  at  nine 
other  standard  pressure  levels  up  to  and  including  100  mil- 
libars (Table  I  (a)).   Those  soundings  were  then  modeled 
into  a  radiative  sounding  model  compatible  in  the  vertical 
with  FNWC  primitive  equation  prediction  levels  (Fig.  2) . 

The  Pacific  Ocean  gridpoint  data  cases  were  taken  from 
the  0000GMT,  16  January,  1974  FNWC  analyses,  while  the 
Atlantic  Ocean  data  cases  were  taken  from  the  1200GMT,  16 
January  1974  analyses.   The  reason  for  the  choice  of  the 
0000GMT  time-analysis  was  to  correspond  to  a  set  of  "sound- 
ings" as  close  as  possible  to  local  solar  noon  over  the  mid- 
Pacific,  whereas  those  of  1200GMT  approximate  a  near-noon 
conditions  over  the  Atlantic  meridian. 

At  each  gridpoint  selected,  the  original  data  listed  the 

five  moisture  values  given  dewpoint  depressions  for  the 

analysis  levels  from  925  mb  to  400  mb .   Since  the  standard 

instrument  level  vapor  pressure  (e  -  )  was  missing  from  the 

ax  r 

original  FNWC  analysis  at  all  gridpoints,  e  .   was  approxi- 

air 

mated  using  the  FNWC  field  of  e  ,  a  computed  value  of  vapor 
pressure  at  about  20  meters  above  MSL  in  the  turbulent 
boundary  layer.   The  FNWC  values  of  T   and  e   over  the  ocean 
had  been  made  available  using  an  operational  planetary  bound- 
ary layer  model  detailed  by  Kaitala  [1974] . 

In  Sections  V  and  VI  it  was  necessary  to  recompute  values 
of  e   (and  of  T  )  by  application  of  the  FNWC  planetary  bound- 
ary  layer  model  to  data  reported  in  the  gridpoint  soundings. 
It  was  found  that  FNWC  values  of  e   as  initially  reported  were 
accurate  enough  to  be  representative  of  the  near-surface 
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Table  1(a).   Example  of  typical  FNWC  sounding  for  gridpoint 
(1,1).   The  humidity  parameters  between  925, . . . , 400mb  are 
dew-point-  depressions,  and  at  the  surface  and  top  of  the 
constant  flux  layer  (9999*)  are  vapor  pressures. 


Pressure 

(mb) 

T(°C) 

surface, 

lOOOmb 

25.8 

925 

18.7 

850 

14.9 

700 

5.9 

500 

-10.7 

400 

-22.4 

300 

-37.3 

250 

-46.2 

200 

-56.2 

150 

-67.4 

100 

-80.7 

9999* 

T  =22.9 

X 

Humidity  Parameters 


ex=24. 

,9 

mb 

1, 

.  5C 

'C 

1, 

.6' 

'C 

4 

.4C 

*C 

4 

.9< 

'C 

2 

,0l 

>C 

*Code  9999  indicates  data  taken  from  the  top  of  the  constant- 
flux  level  (CFL)  of  the  FNWC  initial  data  program. 
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Table  1(b).   Example  of  the  corresponding  radiative  sounding 
with  mixing  ratio  listed  at  odd  k-levels  (Fig.  2).   Addition- 
ally, water  vapor  and  C02  absorber  masses,  cloud  amounts, 
CL(1)  and  CL(2)  where  applicable  in  the  radiative  sounding. 


ssure 

(mb) 

T(°K) 

Mixing  Ratio 

Absorber 

Masses 

Cloud  Amounts 

g/kg 

Water  Vapor 
(gm/cm2) 

C02 
(cm/ cm2) 

CL(1),CL(2) 

1000 

25.8 

15.49 

900 

17.1 

12.76 

CL(2) 

800 

12.1 

2.392 

45.53 

1.000 

700 

5.9 

6.12 

600 

-  1.3 

3.349 

83.53 

CL(1) 

500 

-10.7 

2.31 

.650 

400 

-22.4 

1.37 

3.633 

113.35 

300 

-37.3 

.54 

200 

-56.2 

3.679 

134.00 

100 

-80.7 

.02 

3.680 

138.67 

0 

-80.7 

.00 

3.680 

143.35 
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humidity  values,  e  .  ,  even  though  the  computations  for  e 
made  in  Sections  V  and  VI  employed  other  parameter-values 
for  these  computations. 

In  addition  to  the  parameters  listed  in  Table  I (a,  b) , 
surface  geostrophic  windspeeds  were  also  calculated.   In  the 
procedure  for  calculation  of  the  windspeeds  at  the  gridpoints, 
use  was  made  of  the  100  mb  FNWC  analysis  to  determine  the 
contour  heights  as  D-values.   At  certain  locations,  climatic 
geostrophic  windspeeds  were   used   (Atkinson,  1970).   A 
more  complete  explanation  of  the  use  of  the  geostrophic 
windspeed  calculations  and  its  use  in  connection  with  the 
planetary  boundary  layer  sensible  and  latent  heat  transports 
is  found  in  Section  V,  and  the  windspeed  values  are  listed 
in  that  section. 

To  perform  the  radiative  calculations  indicated  in  Table 
I  (b) ,  data-processing  methods  were  applied  to  convert  both 
the  temperature  and  mixing-ratio  data  at  k-levels  to  obtain 
the  cloud  amounts,  CL  (1)  and  CL  (2)  at  certain  required  k- 
level  boundaries  (Fig.  £)  .   It  should  be  observed  that  all 
soundings  in  this  study  start  at  sea  level  with  the  approxi- 

• 

mation  of  surface  pressure  it  =  1000  mb .   Thus  the  k-levels 
correspond  approximately  to  the  FNWC  a-levels  and  are  taken 
as  P  =  1000.,  900.,  800.,  ....,  200.,  100.,  0.0  respectively. 

A.   INTERPOLATIVE  PROCESSING  TO  k-LEVELS  IN  RADIATIVE  SOUNDINGS 
1 .   Temperatures 

At  each  mandatory  level  of  Table  1(a)  between  the 
1000.,...,  100  mb ,  gridpoint  temperatures  were  listed.   An 
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assumption  was  made  that  at  the  top  of  the  atmosphere,  the 
temperature  was  isothermal  from  100  mb  to  p  =  0.0  mb .   The 
listed  FNWC  sea-surface  temperature  was  taken  as  T10.   At 
the  remaining  k-levels,  air  temperatures  were  derived  from 
either  the  corresponding  FNWC  listed  temperature  or  where 
necessary  by  interpolation  from  listed  values  of  the  original 
FNWC  data.   The  interpolation  scheme  used  was  that  of  the 
three-point  Lagrangian  formula  in  the  vertical. 

Tk  "  To  «pk-piHPk-P2))/((P0-P1)(P0-P1)) 

♦  Tj  ((Pk-P0)(Pk-P2))ACP1-Pe)(P1-P2))       (2-1) 

+  T2  C(Pk-P0)(Pk-P1))/CCP2-P0)(P2-P1)). 

It  is  noted  that  in  Equation  2-1,  k  is  the  required  pressure 

level  needed  for  the  T,  interpolation.   The  level  P   lies 

k      r  o 

above  P,  and  the  pressure  levels  P..  ,  and  P~  lie  below  level 
k.   The  temperatures  T  ,  T- ,  and  T~  are  at  the  same  levels 
as  P  ,  P..  and  P~  respectively. 
2 .   Humidities  At  k-Levels 

Moisture  parameters  from  the  original  gridpoint 
soundings  (Table  I  (a))  were  converted  into  mixing  ratios 
at  each  of  the  original  sounding  levels  (Table  I  (a)). 
Then  utilizing  Equation  2-2  the  near-surface  vapor  pressure 
(taken  as  e  )  was  employed  to  calculate  the  mixing  ratio  at 
k  =  10  as  follows: 

q1000  =  621.97  (eair/1000)  (2-2) 

where  e  .   =  e   is  in  mb  and  qinnn  is  in  gm(kgm)   .   The 
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remaining  mixing  ratios  were  calculated  from  the  original 
FNWC  dewpoint  depression  data  levels  (Table  I  (a))  using 
Equation  2-3  (after  Fleagle  and  Businger,  1963): 

,„>    621.97  exp  (A-B/T)      r  "B^T"TD^  ,       r-  ,. 
q(P)  =  p—^ — L   exP  t  -i^r 1       (2-3) 

where  q(P)  =  mixing  ratio  at  pressure  level  P  and 
A  =  21.656 
B  =  5418. 0°K 

T  =  temperature  (deg.  K)  at  level  P 
T~  =  dewpoint  (deg.  K)  at  level  P 
T-Tn  =  dewpoint  depression  at  level  P. 
The  computed  q-values  were  subsequently  interpolated  to  k- 
levels  using  the  procedure  indicated  in  Equation  2-1.   An 
example  of  the  interpolated  q-values  is  shown  for  the  case 
of  the  radiative  sounding  at  point  (1,  1)  (Table  I  (b)). 
3.   Humidity  Extrapolation  to  p  <  300  mb 

Due  to  the  fact  that  most  radiosonde  humidity  data 
for  p  <  300  mb  is  either  unreliable  or  not  available,  a 
power-law  extrapolation  formula  (Equation  2-4)  similar  to 
one  used  by  Smith  [1966],  was  employed  to  obtain  q-values  at 
k  =  3,  2,  and  1: 

SU21  =  r  .1-  )x  (2-4) 

q5     l  500  J  ^  J 

where  q(P)  =  mixing  ratio  at  level  k  in  g  Kg 

q.  =  mixing  ratio  at  level  P  =  500  mb  in  g  Kg 
P  =  pressure  at  level  k  in  mb . 
Equation  2-4  was  solved  for  the  best-fit  parameter  A,  relating 
the  variables  y  =  log(q(P)/q5)  and  x  =  log(P/500).   Six  q- 
values  from  the  original  sounding  including  P  =  1000,  925, 
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850,  700,  500,  and  400  mb  from  Table  I  (a)  are  used  in  the 
best  fit  solution  (e.g.,  Crow  et   al.,  1960): 


X  = 


iz6 

i-1  yiXi 

i=6 

.2;  x.2 
i=l  i 


(2-5) 


After  the  determination  of  X  from  Equation  2-5,  extrapolation 
of  q-values  from  P  =  400  to  P  =  300,  200,  and  100  mb ,  was 
performed  utilizing  Equation  2-4. 

The  procedure  of  the  previous  paragraph  was  carried 
out  at  each  gridpoint.   As  a  test  of  the  usefulness  of  the 
X-profile  technique  in  determining  q-values  above  300  mb ,  a 
correlation  coefficient  was  statistically  computed  for  each 
X-profile.   The  correlation  coefficient  computed  was  between 
y  and  x  and  had  the  form: 


(2-6) 


As  was  found  in  the  studies  of  Jenks  [1974]  and  Warner  [1974], 
the  correlation  coefficient  was  found  generally  to  lie  in 
the  range  of  .95  to  .99.   Such  high  correlations  indicate 
that  the  X-profile  technique  for  determining  upper  atmosphere 
q-values  has  realistic  estimation  value  for  the  vertical-scale 
involved. 

4 .   Pressure-Scaled  Absorber  Masses 

By  using  the  mixing-ratio  values  computed  at  the 
five  odd  k-levels,  the  layer  pressure-scaled  water  vapor 
absorber  masses  denoted  by  AU  (Fig.  l)  were  calculated  for 
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each  layer  centered  at  the  five  odd  k-levels.   The  general 
equation  for  computing  the  pressure-scaled  water  vapor  mass 
was : 

qkAP     Pk      72 
U(k-l,k+l)  =  JL_  [  m|_r72  (2.?) 

where  the  pressure  ratio  raised  to  the  exponent  0.72  has 
been  employed  (after  Mb'ller  and  Raschke,  1964).   k  is  the 
odd-level  centered  in  each  vertical  layer  concerned  and 
AP  =  200  mb  was  used  for  each  of  the  five  layers.   The  pro- 
cedure was  followed  for  each  layer  (8,  10),  (6,  8),  (4,  6), 
(2,  4) ,  and  (0,  2) . 

The  algorithm  for  computing  the  integrated  water 
vapor  mass  above  the  earth's  surface  was  then  developed  as 
follows : 

U(10)  =  0.0 
U(8)  =  U(10)  +  U(8,  10) 
U(6)  =  U(8)  +  U(6,  8)  / 

U(4)  =  U(6)  +  U(4,  6)  (2-8) 

U(2)  =  U(4)  +  U(2,  4) 
U(l)  =  U(2)  +  0.5U(0,  2) 
U(0)  =  U(2)  +  U(0,  2). 

In  the  computation  of  the  carbon-dioxide  scaled 
absorber  mass,  a  computational  routine  similar  to  the  one 
used  for  water  vapor  was  employed.   Because  it  is  customary 
to  state  the  C0?  absorber  masses  in  terms  of  the  N.  T.  P. 
volume  of  the  gas  in  cms  over  a  cm2  of  the  earth's  surface, 
the  C02  absorber  mass  (AC)  centered  on  an  odd  k-level  becomes 
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(2-9) 
ACk  =  C(k-l,k+l]  =  3.14  x  10"*  (  |£-  )(  T5l|-j5  )0-72. 

Here  k  is  the  odd-level  indicator,  3.14  x  10    is  the  con- 
stant mixing  ratio  of  CO-  by  volume,  while  AP/gp   is  the  N. 
T.  P.  thickness  of  any  layer  of  pressure- increment  AP  =  200 
mb.   p   is  the  constant  density  in  a  homogeneous  atmosphere 
of  depth  H,  in  which 

TrSurfl  =  273-16  de§-  K  and  ps  =  1013-25  mb- 

Finally,  Equation  2-9  is  shown  to  be  represented  in  the  form 

p 

C(k-l,k+l)    -    3.14   x    10"*CH)(   101^25    )(   T0l3725    >  *  " 

(2-10) 
RT 

H   =   — -  =    7.9954   x    105    cm. 
g 

Equation  2-10  was  then  utilized  for  each  odd  level  k  =  9,  7, 
5,  3,  1. 

As  was  done  in  the  water-vapor  case,  an  algorithm 
similar  to  Equation  2-8  is  developed  to  find  the  integrated 
carbon  dioxide  mass  above  the  earth's  surface: 

C(10)  =  0.0 

C(8)  =  C(10)  +  C(8,  10) 

C(6)  =  C(8)  +  C(6,  8) 

C(4)  =  C(6)  +  C(4,  6)  (2-11) 

C(2)  =  C(4)  +  C(2,  4) 

C(l)  =  C(2)  +  0.5C(0,  2) 

C(0)  =  C(2)  +  C(0,  2)  . 
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B.   CLOUD  PARAMETERIZATION 

The  relative  humidities  and  thus  the  saturation  vapor 
pressure  at  levels  k  =  5  and  k  =  9  are  used  in  the  calcula- 
tions of  the  fractional  cloud  cover  amounts  in  layers  (4,  6) 
and  (8,  9).   The  saturation  vapor  pressures  at  each  level 
were  computed  using  the  following  equations  (cf . ,  Equation 

2-3)  : 

621.97  exp  [  A  -  ^_  ] 

qs(5)  =  ^ S- (2-12  (a)) 

621.97  exp  [  A  -  ~  ] 
qs(9)  =  900 " C2"12  (b» 

where  A  =  21.656  and  B  =  5418.0  deg.  K. 

The  relative  humidities  were  then  calculated  from: 

R.H.(5)  =  ^|iy  (2-13  (a)) 

R.H.(9)  =  3I9L-  .  (2-13  (b)) 

Using  these  relative  humidity  values,  the  fractional  cloud 
amounts  for  each  layer  were  parameterized  employing  equations 
similar  to  Smagorinsky  [1960] : 

CL(1)  =  2.0  (RH(5))  -  0.7  (2-14  (a)) 

CL(2)  =  3.33  (RH(9))  -  2.0.  (2-14  (b)) 

Equations  2-14  have  been  employed  in  the  radiative  studies  of 
Jenks  [1974]  for  a  spring  day  in  1973  and  by  Warner  [1974]  for 
a  case  in  mid-October,  1973.   In  the  current  study,  Equations 
2-14  are  tested  in  a  winter  sample  of  FNWC  data  for  16 
January  1974. 
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Smagorinsky ' s  parameterization  of  CL(1)  and  CL(2)  per- 
mits cloudiness  fractions  of  1.30  and  1.333  respectively  with 
RH  =  1.0.   Values  of  CL  of  this  magnitude  were  allowed  orig- 
inally by  Smagorinsky  to  suggest  the  type  of  supersaturation 
which  accompanies  precipitation.   However,  in  all  of  the 
radiative  models  in  use  (that  of  FNWC ,  UCLA,  NCAR,  etc.), 
it  was  considered  correct  procedure  to  limit  CL  <  1.0  in 
both  layers.   In  the  radiative  calculations  of  this  study  an 
arbitrary  reduction  by  one  third  of  that  computed  by  the 
Smagorinsky  formulations  (2-14)  was  made  so  that  neither 
fractional  cloud-cover  could  exceed  unity.   For  radiative 
calculations,  the  degree  of  supersaturation  in  clouds  was 
considered  not  of  importance. 

Thus  CL(1)  and  CL(2)  computed  by  (2-l4as  b)  are 
replaced  respectively  by  their  comparison  cloud-covers, 
CL(1)'  and  CL(2)',  defined  by 

CL(1)  '  =  2/3  CL(1) 

(2-15) 
CL(2) •  =  2/3  CL(2) 

computed  at  each  gridpoint.   The  objective  for  testing  the 
reduced  cloud-covers  (2-15)  was  to  determine  the  relative 
sense  of  the  displacement  of  the  radiative  balance  of  the 
earth  atmosphere  system.   Such  radiative  tests  were  especially 
crucial  since  evidence  from  Warner's  [1974]  tests  made  with 
unmodified  (CL(1),  CL(2))  indicated  that  the  planetary  albedo 
was  larger  than  observations  from  NIMBUS  III  (Raschke,  et   al. 
1973).   The  one-third  reduction  in  the  Smagorinsky  calculated 
amounts  used  in  this  study  is  an  initial  effort  to  tune  cloud 
amounts  to  give  radiational  results  in  closer  agreement  with 
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the  latest  satellite  results  such  as  are  found  in  Raschke 
et   at.     [1973].   Also,  the  fractional  cloud  amounts  were 
considered  functions  of  large-scale  effects  only.   Small 
scale  convective  activity,  seasonal  conditions  (other  than 
the  soundings  themselves) ,  and  latitudinal  effects  were  not 
considered  in  specifying  this  reduction  factor.   The  purpose 
of  the  reduced  cloud-parameterization  modification  was  for 
estimating  radiational  effects  only. 

C.   WINDSPEED  USED  IN  TURBULENT  HEAT  TRANSFERS 

For  all  gridpoints  located  between  5°N  and  5°S,  gradient 
windspeeds  were  taken  from  the  climatological  data  for  the 
month  of  January  (Atkinson,  1970).   Gradient  isotach  values 
were  utilized  at  the  tropical  gridpoints  noted  due  to  the 
small  Coriolis  parameter  f  in  Equation  (2- 16) (below)  in  the 
latitudes  between  5°N  and  5°S. 

At  all  other  gridpoints,  Equation  (2-16)  was  utilized  to 
compute  geostrophic  windspeed  given  by 

IV  1  -  f-l^fL  (2-16) 

here  g  =  9.80665m  sec 

r      +  1.86603 

m  =  map  factor  =  r — -  c . — , 
^  1  +  Sine}) 

<J>  =  latitude 

d  =  381.  x  103m  =  nominal  FNWC  grid-dimension 

f  =  <  86TI00  )  Sin*  t2"17) 

I  Az|  =  /(z4  -  z2r  +  Cz3  -  zi)2  .  (2-18) 
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In  (2-16)  and  (2-18),  Az  is  the  contour-difference  cen- 
tered in  general  on  the  gridpoint.  At  the  map  edge  grid- 
points  of  meridians  1,  2,  3,  4  it  was  possible  to  compute 
only  a  forward  differenced  version  of  Equation  2-18  (cf . , 
Fig.  3).  Thus  at  edge-points  "0",  Equation  2-19  was  used 
(below).  Note,  the  grid  distance  d/m  in  Equation  2-19  is 
one-half  of  that  used  in  Equation  (2-16). 

|V  I  =  |  ,hz|  (2-19) 

1  g1    f  (d/m)  v     J 

with 

Az  =  /(z4  -  Zo)2  +  (z3  -  z0)2  .  (2-20) 

All  geostrophic  windspeeds  were  then  converted  to  surface 
windspeeds  using  the  following  empiricism  due  to  Langlois 
and  Kwok,  [1969]  . 

V    =  .8  |V  I  +  2.2  (2-21) 

s  I   g.l  K  J 

where  V     =  surface  windspeed  (mps) 

V   =  1000  mb  geostrophic  windspeed. 

Note  that  in  the  interval  5  to  10  degrees  of  latitude  in 
either  hemisphere,  the  value  of  the  Coriolis  parameter  was 
arbitrarily  set  at  f  =  .25  x  10   sec   ,  following  Kaitala 
[1974],  while  for  |<j>|  >  10°,  the  actual  value  of  f  for  this 
particular  latitude  was  used. 
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Figure  3.   Section  of  FNWC  polar  stereographic  meridian  il- 
lustrating the  method  of  obtaining  contour  gradients  in  the 
vicinity  of  a  gridpoint.   The  values  of  contour-heights 
correspond  to  the  1000  mb  D-values  in  the  neighborhood  of 
the  gridpoint. 
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III.   TERRESTRIAL  RADIATION 

A.   THEORETICAL  AND  EMPIRICAL  BASIS 

Empirical  formulas  for  computing  emissivities  for  ter- 
restrial flux  calculations  in  the  atmosphere  were  developed 
by  Sasamori  [1968]  in  association  with  the  operation  of  the 
NCAR  General  Circulation  Model.   These  empirical  formulas 
for  water  vapor  and  C0~  emissivities  were  matched  to  the 
theoretical  values  used  in  the  Radiation  Chart  developed  by 
Yamamoto  [1952].   This  chart  was  proved  by  Jenks  [1974]  to 
be  accurate  enough  to  simulate  the  empirical  methods  of 
Sasmori  and  is  used  here  as  a  schematic  guide  for  the  inte- 
gration of  the  radiative  transfers  through  the  various  layers 
in  the  radiative  soundings  described  in  Section  II. 

For  use  in  the  FNWC  heating  package,  the  essential  long- 
wave flux  formulas  required  are  the  following: 

F   *  =  net  IR  flux  at  earth,  k  =  10 
F,*  =  net  IR  flux  at  level,  k  =  6 
F~*  =  net  IR  flux  at  level,  k  =  2 
F610  =  net  IR  flux  divergence  in  the  layer  (6,  10) 
F26  =  net  IR  flux  divergence  in  the  layer  (2,  6). 
In  order  to  get  the  IR  flux-divergences  (i.e.,  cooling 
rates)  in  the  layers  (6,  10)  and  (2,  6),  the  differences 
Fg*  -  F^q*  and  F  *  -  Ffi*  must  be  computed.   The  detailed 
scheme  for  making  such  computations  for  various  combinations 
of  cloud  cover  CL(1),  CL(2)  is  similar  to  that  used  by  Warner 
[1974]  and  is  explained  in  the  following  subsections  B,  C,  D, 
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and  E  along  with  the  diagrams  Figures  (4a,b,c,d),  (5a,b,csd) 
and   (6a,b,c,d).  A  complete  listing  of  net  IR  fluxes  and 
cloud  amounts  was  computed  at  each  gridpoint.   An  example 
of  this  listing  is  shown  for  "full-cloud"  amounts,  Equation 
2-14,  in  Table  II.   One  must  further  have  a  physically  sound 

representation  of  the  emissivity  e    as  a  function  of  both 

r  J      wc 

water-vapor  and  C02  absorber  masses  along  the  sounding,  that 
is,  along  the  path  of  integration  in  order  to  make  IR  net- 
flux  calculations.   For  further  explanation  of  the  emissivity 
formulas  (after  Sasamori  1968)  used  in  the  quadrature  scheme, 
refer  to  Appendix  A. 

B.   NET  FLUX  F1Q* 

As  depicted  in  Table  I  (b)  ,  radiative  soundings  have  been 
computed  as  the  combination  of  parameters  U(k,  10),  C(k,  10), 
T, ,  for  each  required  level  k.   Also,  listed  at  each  grid- 
point  are  the  cloud-cover  parameters  CL(1)  and  CL(2)  com- 
puted from  (2-14).   The  grid  area  may  then  be  thought  of 
being  composed  of  areal  fractions  (or  weights)  W(0,  0), 
W(l,  1),  W(l,  0),  and  W(0,  1)  which  are  defined  by  Equations 
3-1, >3-4  as  follows: 

(a)  totally  clear  skies,  case  (0,  0) 

W(0,  0)  =  (1-CL(1))(1-CL(2))  (3-1) 

(b)  overcast  in  both  layers,  case  (1,  1) 

W(l,  1)  =  (CL(1))(CL(2))  (3-2) 

(c)  overcast  in  upper  layer  only,  (case  1,  0) 

W(l,  0)  =  (CL(1))(1-CL(2))  (3-3) 

(d)  overcast  in  lower  layer  only,  (case  0,  1) 

W(0,  1)  =  (1-CL(1))(CL(2)).  (3-4) 
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By  multiplying  the  appropriate  weighting  factors  by 
the  respective  net  flux  computations  F..  *(  0,  0),..., 
F-.Q*(0,  1)  the  composite  net  fluxes  at  level  k  =  10  for 
the  combined  cloud  (CL(1),  CL(2))  coverage  over  the  grid 
area  are  then  obtained.   Such  composite  net  flux  values 
have  the  following  form  using  level  k  =  10  as  an  example: 

F1Q*  =  W(0,  0)F1Q*(0,  0)  +  W(l,  1)F10*(1,  1) 

+  W(l,  0)F10*(1,  0)  +  W(0,  1)F10*(0,  1). 

1.   Net  Flux  with  Clear  Skies,  F10*(0>  0) 

In  Fig.  4  (a)  the  hatched  area  represents  downward 
flux  F10i  at  k  =  10  or  in  integral  form: 


(3-5) 


B10 

F   4-  = 
WC     B^O 


J  e    dB  ,   Bk  =   aTk 


Bv  is  the  Stephen-Boltzman  blackbody  flux  at  T,  .   e    is  the 

JC  K      WC 

combined  water-vapor  and  CC"   emissivity  along  the  sounding, 

and  is  related  to  absorber-masses  as  detailed  in  Appendix  A. 

The  unhatched  area  therefore  depicts  the  net  flux  at  k  =  10 

with' clear  skies  given  in  the  equation 

^10 
F10*(°>  °)  •  A  U  "  £wc)dB'  B10  =  Stnio"-     (3"6) 

D-U 

2 .   Net  Flux  with  Overcast  Clouds  in  Both  Layers,  F-.  *(!,!) 
Figure  4  (b)  illustrates  this  case.   By  using  the 
trapezoidal  summation  techniques,  F10I  at  k  =  10  is  depicted 
once  again  by  the  hatched  area  representative  of  the  case  of 
an  overcast  with  base  at  k  =  9 .   The  unhatched  area  depicts 
the  integral 


^5 


F   * 


plO 
(1,  1)  -   J     (1  "  ewc)dB  (3-7) 

B=0 
with  overcast  clouds  at  both  upper  and  lower  levels  (Fig.  4 

(b)). 

3 .  Net  Flux  with  Overcast  Clouds  in  Upper  Layer  Only, 

EIqSL  o) 

Here  in  Fig.  4  (c) ,  the  hatched  area  depicts  the 
total  F1(H  emitted  down  from  the  cloud  base  at  k  =  6,  plus 
an  increment  in  the  upper  right  corner  of  Fig.  4  (c) ,  which 
represents  the  contribution  to  F4-  from  the  water  vapor  and 
CO?  in  the  atmospheric  layer  from  k  =  6  to  k  =  10.   The  un- 
hatched  area  depicts  the  integral 

plO 
F10*C1,  0)  =  J  (1  -  ewc)dB.  (3-8) 

D— (J 

4.  Net  Flux  with  Overcast  in  Lower  Layer  Only  F-  *(0,  1) 
This  case  is  the  same  as  depicted  for  case  F..  *(1,  1) 

in  Fig.  4  (d)  and  is  explained  in  Section  III.(B).2.   e    is 

the  temperature-dependent  emissivity  for  the  completion  of 

the  integration  from  level  k  =  1  to  the  origin  (B  =  0 ,  e  = 

1.0)  of  the  Yamamoto  Radiation  Chart,  applied  for  T  <  210 

deg.  K.   e    is  the  temperature-independent  emissivity  applied 

where  T  >  210  deg.  K  in  the  level  from  k  =  1  to  k  =  10.   For 

the  specific  formulations  of  e    and  e    after  Sasamori  [1968] 
r  wc      wc 

see  Appendix  A. 

5.  The  Composite  F*  (CL(1)  CL(2))  Equation 

The  principle  involved  in  compositing  net  fluxes  has 
already  been  expressed  by  Equation  3-5.  The  following  equa- 
tion can  be  thought  of  as  the  general  form  of  the  quadrature 


k6 


scheme  used  to  evaluate  the  integrals  for  computing  F..  *  for 
any  combination  of  (CL(1),  CL(2)).   It  is  also  used  in  com- 
puting the  integral  for  the  composite  cloud-cover  case  F. 
corresponding  to  (3-5): 


* 

0  * 


F10*(CL1,CL2)  =  [l-CL(2)]{(B10-B6)-.5[cwc(8,10)(B10-B8) 

+  (e   (8,10)+e   (6,10))  (B0-B,)]} 

+  (l-CL(2))(l-CL(l)){B6-.5[(£wc(6,10) 

+  cwc(4,10))(B6-B4)+(£wc(4,10)         (3-9) 

+  ewc(2,10))(B4-B2)+(ewc(2,10) 

+  ewc(l,0))(B2-B1)+ew(;((0,10),T1)*B1]} 

+  CL(2)((B10-Bg)[l-.5ewc(9,10)]}. 

By  setting  CL(1)  =  CL(2)  =  0.0,  Equation  3-9  reduces 
to  the  form  needed  in  calculating  F10*(0,  0).   Similarly,  one 
can  obtain  the  equation  for  computing  F-  *(1,  1)  as  depicted 
on  the  chart  in  Fig.  4  (b)  by  setting  CL(1)  =  CL(2)  =  1.0. 
Finally,  the  other  two  subcases  F..  *(1,  0)  and  F-fi*(0,  1)) 
can  be  formulated  by  setting  CL(1)  =  1.0,  CL(2)  =  0.0,  and 
CL(1)  =  0.0,  CL(2)  =  1.0,  respectively  in  Equation  3-9. 
Similar  procedures  were  followed  for  computing  net  IR  flux 
at  levels  k  =  6  and  k  =  2,  as  described  in  the  Subsections 
C  and  D  (Figs.  5  (a),  (b) ,  (c) ,  (d)  and  6  (a),  (b) ,  (c) ,  (d) , 
respectively. 
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C.   NET  FLUX  F  * 
o 

1 .  The  Composite  F  *  Equation 

The  following  expression,  which  represents  the  F  * 
composite  net  flux,  may  be  thought  of  as  the  weighted  sum  of 
the  unhatched  areas  of  Figs.  5  (a),  (b) ,  (c) ,  and  (d) ,  with 
the  appropriate  weight  factors  W(0,  0),  W(l,  1),  W(l,  0), 
and  W(0,  1)  of  Equations  3-1,  2,  3,  and  4  applied.   The  re- 
sulting expression  for  F,*(CL..,  CL?)  becomes  (after  Warner, 
1974) 

F6*  =  [l-CL(l)]{B8-.5[ewc(6,8)(B8-B6)+ewc(4,6)(B6-B4) 

+  (e   (4,6)+e   (2 ,6) ) (B, -B0) + (e   (2,6)+e   (1 , 6) ) (B~ -B- ) 

^  wcv  '  **   wcv  '  JJK    4  2J     K   wc^  '  }      wc  *•  '  JJK    2      1J 

+  ewc((0,6),T1)*B1]}+(l-CL(l))(l-CL(2)){(B10-B8)[l. 
-  .5(ewc(6,8)+cwc(6,10))]}+CL(l){(B8-B6)*        (3-10) 
[l-.5ewc(6,8)]}+CL(l)(l-CL(2)){(B10-B8)* 
[l-.5(ewc(6,8)+ewc(6,10))]}. 

2.  Net  Flux  With  Clear  Skies  F6*(0,  0) 


As  was  described  in  Section  III.B.5  concerning  the 
use  of  the  composite  case  (Equation  3-9),  Equation  3-10  can 
be  reduced  to  give  expressions  for  F,.*  for  the  various  cloud- 
cover  cases  (0,  0),  (1,  1),  (1,  0),  and  (0,  1).   The  resulting 
expressions  are  depicted  as  the  unhatched  areas  in  Figs.  5  (a), 
(b) ,  (c) ,  and  (d) .   For  the  clear  sky  case  of  Ffi*,  Equation 
5-10  would  be  employed  with  CL(1)  =  CL(2)  =  0.0.   This  is 
depicted  by  the  unhatched  area  in  Fig.  5  (a) . 
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3 .  Net  Flux  With  Overcast  Clouds  in  Both  Layers  F  *(1,  1) 
For  this  case,  Equation  3-10  would  be  used  with  CL(1)  = 

CL(2)  =  1.0.   This  is  depicted  by  the  unhatched  area  in  Fig. 
5  (b). 

4.  Net  Flux  With  Overcast  Clouds  in  Upper  Layer  Only 

The  expression  needed  in  this  case,  and  which  is 
depicted  by  the  unhatched  area  in  Fig.  5  (c) ,  is  formulated 
by  using  Equation  3-10  and  setting  CL(1)  =  1.0  and  CL(2)  = 
0.0. 

5 .  Net  Flux  With  Overcast  Clouds  in  Lower  Layer  Only 
F6*(0,  1) 

The  expression  needed  in  this  case,  representable  by 
the  unhatched  area  in  Fig.  5  (d) ,  is  formulated  by  using 
Equation  3-10  and  setting  CL(1)  =  0.0  and  CL(2)  =  1.0. 

D.   NET  FLUX  F 2* 

A  procedure  analogous  to  the  ones  used  in  the  computa- 
tions of  F-.  *  and  F,-*  net  fluxes  can  be  used  in  calculation 
1U       o 

the  net  flux  at  level  k  =  2.   F7*  is  represented  pictorially 
by  the  weighted  sum  of  the  unhatched  areas  of  Figs.  6  (a), 
(b) ,  (c) ,  and  (d) .   Once  these  areas  have  been  computed  with 
the  proper  weighting  factors  applied,  as  was  done  similarly 
in  Equations  3-9  and  3-10,  the  following  general  expression 
for  the  composite  result  for  F~*  results 

F2*  =  [l-CL(l)]{B8-.5[ewc(2,4)(B4-B2)+(ewc(2,4) 

+  S,c(2,6))(B6-B4)+(£wc(2,6)+ewc(2,8))* 

(B8"B6)+ewc(1'2)(B2"Bl)+£wc((0'2)'Tl)*Bl)]}      (3'11} 
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+  (l-CL(l))(l-CL(2))*{(B10-B8)[l-.5(ewc(2,8)       (3-11) 
+  ewc(2,10))]}+CL(l){B4-.5[£wc(2,4)(B4-B2) 
+  ewc(l,2)(B2-B1)+ewc((0,2),T1)*B1]} 

Equation  3-11  can  also  be  derived  in  the  computation  of  the 
net  fluxes  for  the  general  cloud-cover  cases  once  the  proper 
weight  factors  for  the  simplified  net  fluxes  are  introduced. 
These  simplified  F-*  cases  correspond  to  the  (0,  0),  (1,  1), 
(1,  0),  (0,  1)  cloud-overcast/clear  combinations  and  are 
represented  by  the  notation  F  *(0,  0),  F  *(1,  1),  F  *(1,  0), 
F  *(0,  1),  respectively.  The  corresponding  chart  depictions 
are  displayed  in  Fig.  6  (a),  (b) ,  (c) ,  (d) ,  respectively. 

E.   APPLICATION  TO  HEAT  BALANCE  COMPUTATIONS 

In  order  to  compute  the  heat  balance  of  the  atmosphere, 
layers  (2,  6)  and  (6,  10)  are  treated  as  subject  to  solar 
heating  (Section  IV) .   In  order  to  complete  the  atmospheric 
radiational  balance,  the  long-wave  mean  cooling  rate  in  the 
FNWC  model  was  computed  as  follows  (after  Kaitala  [1974]): 

F26  =  F2*  -  F6*  (3-12(a)) 

F610  =  F6*  -  F1Q*  (3-12(b)) 

where  the  combination  26  and  610  indicates  the  layer  boundaries. 
Equations  3-12 (a)  and  (b)  represent  the  cooling  effects  by  ter- 
restrial radiation  in  the  indicated  layers,  which  is  applicable 
at  the  time  of  the  analysis.   The  radiative  balance  at  the  sur- 
face is  computed  using  the  absorbed  solar  insolation  (Section 
IV.B-C)  and  the  net  flux  F-  *  from  Equation  (3-9).   A  complete 
heat  balance  at  the  surface  is  computed  by  subtraction  of  the 
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evaporative  and  sensible  heat  transport  (Section  V)  from  the 
surface  radiative  "balance."  All  these  balances  are  dis- 
played in  cross  sectional  display  forms  in  Section  VI.   The 
IR-cooling  analysis  was  carried  out  at  each  gridpoint  and 
the  cooling  rates  were  stored  by  gridpoint  along  each  meridian 
The  zonally-averaged  atmospheric  balances  are  detailed  in 
Sections  VII  and  VIII  for  16  January  1974. 

F.   STATISTICAL  RESULTS  AND  COMPARISONS 

A  statistical  test  of  the  F-j0*  numerical  results  (Equa- 
tion 3-9)  was  conducted  by  utilizing  a  linear  regression 
program  for  regressing  the  93  values  of  F..  *  on  the  cor- 
responding sample  values  of  Bin  and  B.^/i"  as  predictors.   The 
program  used  was  based  on  the  BMD02R  in  the  Biomedical  set 
of  programs  (after  DIXON,  1973) .   The  resulting  form  was 

F10*  =  B10  (a  +  b/^  C3"13) 


where  Blf)  =  Stephen-Boltzman  blackbody  radiation  at  the  sur- 
face with  temperature  T1Q.   B..^  =  St  *  ^■\(]'* 

e  =  surface  layer  vapor  pressure  (mb)  . 
a  and  b  in  Equation  3-13  are  the  desired  regression  coeffi- 
cients and  are  dependent  on  B1f)  and  B1f)/e~  (after  Brunt,  1932) 
The  predictand  used  for  the  regression  was  Fin*>  which  in- 
cluded the  downward  contribution  of  CO-  as  well  as  that  of 
water  vapor  in  the  computation  of  net  flux  at  level  10.   The 
statistical  equation  found  in  this  study  is  Equation  3-14 

where  R   is  the  multiple  correlation  coefficient, 
m  r 

F   *  =  .3564Bin  -  .0224Bin/e 
±u  iu         iu  (3-14) 

R   =  .9955. 
m 
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The  high  value  of  the  R   shows  that  Birv  and  Bin/e~  are  good 
&  m  10      10       & 

predictors  of  F-fl*  and  therefore  of  F,  obtained  from 
F-Q*  =  B.q  -  F,.   Introduction  of  F..  *  into  3-14  permits 
solution  of  the  latter  for  the  downward  flux  at  the  surface. 
F,  is  given  by 

Fd  =  B1(){(l-a)-b/e"} 
or  (3-15) 

Fd  =  Bio^-6436  +  -0224/e) . 
Equation  3-15  is  to  be  compared  with  a  similar  result  after 
Swinbank  [1963] 

Fd  =  Bl0^'64  +  -O37^)-  (3-16) 

In  summary,  the  high  value  of  the  correlation  coefficient  of 
Equation  3-14  along  with  a  strong  similarity  between  Equations 
3-15  and  3-16  seems  to  indicate  that  Equation  3-14  may  be 
regarded  as  an  oceanic  version  of  the  Brunt  downward-flux 
equation.   Also,  it  is  noted  that  F_,  depends  most  strongly 
upon  B1f)  and  only  slightly  on  (Bin/e~).   This  was  to  be  ex- 
pected over  the  oceanic  regime  where  nearly  constant  relative 
humidities  exist  over  large,  expansive  areas. 

A  second  statistical  test  was  based  upon  a  possible  re- 
lationship between  cloudy  sky  and  clear  sky  cases  of  F-  *. 
In  this  test,  the  ratio  F, n* (CL) /F1 n* (0 ,  0)  was  used  as  the 
predictand.   The  predictor  was  an  approximation  to  the  com- 
posite total  opaque  cloud  at  each  gridpoint,  defined  as: 

CL  =  CL(1)  +  CL(2)(1  -  CL(1)).  (3-17) 

The  regression  was  run  twice,  first  utilizing  F.)f*(CL)  com- 
puted using  the  amounts  of  CL(1)  and  CL(2)  as  calculated  from 
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the  Smagorinsky  formulas  (Equations  2-l4(a),  (b))  and  second, 
utilizing  F..  *(CL)  computed  from  the  two-thirds  amounts  of 
CL(1)  and  CL(2).   The  resulting  regression  formulas  based 
upon  a  best-fit  to  the  93  soundings  used  in  this  study  were 
listed  in  (3-18,  3-19) 


FUL 


L  CL:   F1Q*(CL)  =  F10*(0,  0)(1-.75CL),  R  =  .9916   (3-18) 


2/3  CL:   F1Q*(CL)  =  F10*(0,  0)(1-.73CL),  R  =  .9936.  (3-19) 

A  comparison  equation  after  Sverdrup  [1942]  is  included  as 
(3-20) 

F1()*(CL)  =  F10*(0,  0)(1-.83CL).  (3-20) 

The  high  correlation  coefficients  found  with  Equation  3-18  and 
3-19  along  with  the  similarity  to  the  Equation  3-20  (after 
Sverdrup,  1942),  indicates  the  general  capability  of  the  IR- 
radiative  model  employed  in  this  study  to  account  for  the 
effective  net  radiation  at  the  surface,  F.  *  (CL. , CL~) ,  both 
in  the  completely  clear  and  cloudy-sky  cases. 

The  mean  values  of  the  parameters  in  Equations  3-18  and 
3-19  over  the  93  gridpoint  soundings  were  as  follows: 
FULL-CLOUD  STATISTICS  REDUCED-CLOUD  STATISTICS 

F10*(CL)  =  .0749  ly   min"1  F10*(CL)  =  .0950  ly   min"1 

F10*(0,  0)  =  .1518  £y  min"1         F1Q*(0,  0)  =  .1518  ly   min"1 

CL  =  .6951  CL  =  .5190 

In  the  above  table,  superior-bar  symbols  denote  sample  means. 
These  results  indicate  that  the  surface  net  flux  is  decreased 
by  the  respective  ratios     (. 0769/ . 1518)  =  .507  with  CL  = 
.6951  and  (. 0568/ . 1518)  =  .3740  with  CL  =  .5190.   These  results 
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indicate  that  the  downward  flux  Finl(CL)  is,  on  the  average, 
50.7  percent  larger  than  the  clear-sky  case  utilizing  the 
cloud-parameterization  model  (Equations  2-14(a),  (b))  in 
this  study  with  full  cloud  amounts.   On  the  other  F1f)4-(CL) 
for  the  two-thirds  cloud-model,  the  downward  flux  is  increased 
by  only  37.4  percent  relative  to  the  clear-sky  case. 

G.   COMPARISONS  WITH  CLIMATOLOGICAL  RESULTS  OF  F2 

Comparison  was  made  of  F?  by  the  model-computations  with 
satellite  measurements  of  total  long-wave  flux  to  space  after 
Raschke  et   at.     [1973]  for  the  NIMBUS  III  period  21  January  - 
3  February  1970.   F-  is  used  instead  of  F  *  since  the  satel- 
lite does  not  measure  the  downward  flux  at  k  =  2 .   Fz is 
computed  using  Equation  3-11,  which  was  also  used  for  F  * 
but  omitting  those  emissivity- terms  involving  levels  above 
k  =  2,  that  is  e   (2.  1)  and  z      (2,  0),  from  Equation  3-11. 

i  wc  v.   >    J  wc  v.   ,    j   ,  n 

These  terms  are  associated  with  downward  flux  at  k  =  2 .   The 
results  for  16  January  1974  were  computed  both  for  the  full- 
cloud  and  the  two-thirds  cloud-amount  cases,  based  upon  all 
gridpoint  soundings.   These  F~  -  values  were  then  interpolated 
to  whole  multiple  of  5°  latitude  along  each  of  the  four  merid- 
ians, and  the  results  lat itudinally  averaged  to  get  a  "zonally- 
averaged"  value  valid  for  each  5°  latitude  band  used  in  this 
study. 

Table  III  (a)  shows  the  zonally-averaged  F~  -  values  com- 
pared with  those  after  Raschke  et   at.       For  consistency  in 
comparison  of  the  F~  results  deduced  here  with  those  reported 
by  Raschke,  the  latter's  results  were  obtained  by  interpola- 
tion from  charts  using  the  same  oceanic  meridians  as  those 
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Lat 


Radiative  Statistics 

F2        F2       F2 
(FULL-CL)   (2/3-CL)  (Raschke) 


Cloud  Statistics 

FULL-Cloud   2/3-Cloud 
CL  AVG.     CL  AVG. 


20s 

.311 

.325 

.370 

.932 

.758 

15 

.321 

.334 

.378 

.927 

.733 

10 

.344 

.354 

.377 

.867 

.634 

5 

.354 

.363 

.347 

.746 

.536 

0 

.361 

.369 

.351 

.790 

.556 

5 

.371 

.378 

.368 

.766 

.527 

10 

.382 

.387 

.388 

.557 

.376 

15 

.386 

.393 

.394 

.532 

.342 

20 

.390 

.388 

.391 

.443 

.279 

25 

.361 

.362 

.371 

.463 

.320 

30 

.354 

.328 

.344 

.559 

.411 

35 

.336 

.342 

.319 

.516 

.382 

40 

.315 

.319 

.304 

.628 

.529 

45 

.278 

.270 

.292 

.833 

.740 

50 

.264 

.280 

.270 

.911 

.759 

55 

.217 

.233 

.255 

.931 

.749 

60 

.209 

.226 

.258 

.764 

.528 

65N 

.235 

.235 

.230 

0.0 

0.0 

Wt. 
Avg. 

.333 

.338 

.344 

.688 

.514 

Table  III (a).   Comparison  of  zonally-averaged  longwave  flux 
to  space,  F2,  as  found  by  this  study  for  both  full-CL  and 
2/3-CL  cloud  cases  for  16  January  1974,  and  by  Raschke  et  al. 
(1973)  based  upon  NIMBUS  III  measurements.   Also  included  are 
composite  cloud  amount  fractions  used  in  the  two  cloud 


models  of  this  study.   Flux  values  in  ly  min 
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55 
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used  in  this  study.   For  further  insight  into  the  compari- 
sons, the  cloud  statistics,  CL  and  2/3-CL  values  are  listed 
at  each  latitude,  where  CL  is  given  by  Equation  3-17.   In 
the  bottom  line  of  each  column  in  Table  III  (a)  is  listed 
the  cosine-weighted  mean  of  each  set  of  the  column  values 
over  the  listed  latitude  range.   The  resulting  oceanic-means 
of  F2  obtained  by  the  various  techniques  are  as  follows: 

F2(Full  cloud  amount  this  study)  =  .333  £y  min 
F2(2/3  cloud  amount  this  study)  =  .338  £y  min 
F2 (after  Raschke)  =  .344  Zy   min  *. 

Note  that  the  Raschke  results  do  not  correspond  to  speci- 
fic known  CL  -  values.   The  only  conclusion  that  can  be  drawn 
from  Table  III  (a)  is  that  both  cloud-amount  systems  give  F~ 
results  which  are  reasonably  close  to  the  values  reported  by 
Raschke.   This  is  especially  true  in  the  subtropics,  15°-20°N, 
where  peak  values  of  F~  are  found,  and  in  high  latitudes  where 
minimum  values  of  F~  (after  Raschke,  1973),  are  found.   By 
these  results  alone,  neither  the  full-cloud  model  nor  the 
2/3  CL  can  be  selected  as  preferable.   The  limitations  of 
the  comparisons  made  here  are  obvious,  when  it  is  realized 
that  at  latitudes  between  20°S-5°S  and  between  60°N-65°N  in 
Table  III  (a),  fewer  than  4  meridional  lines  enter  into  the 
zonal  averages  listed.   For  all  other  latitudinal  average 
values,  four  meridional  lines  were  used  in  this  averaging. 
These  same  limitations  will  apply  to  Table  III  (b)  below. 
Note  finally  that  the  model  results  F2  (full)  and  F2  (2/3  CL) 
are  based  upon  a  single-day's  observations,  16  January  1974. 
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F10* 

Fi  o* 

Fi  o* 

Lat 

(FULL  CL) 

(2/3  CL) 

(Budyko) 

20S 

.054 

.078 

.091 

15 

.050 

.079 

.092 

10 

.055 

.086 

.083 

5 

.063 

.090 

.085 

0 

.052 

.082 

.084 

5 

.056 

.086 

.080 

10 

.085 

.107 

.086 

15 

.094 

.115 

.094 

20 

.102 

.120 

.102 

25 

.104 

.120 

.098 

30 

.096 

.111 

.099 

35 

.107 

.122 

.102 

40 

.096 

.108 

.109 

45 

.056 

.067 

.098 

50 

.054 

.069 

.099 

55 

.052 

.068 

.102 

60 

.085 

.106 

.104 

65N 

.095 

.095 

.095 

Wt.  Avg. 

.075 

.096 

.094 

Table  111(b).   Comparison  of  zonally-averaged  net  fluxes  at 
the  surface  F10*  as  found  by  this  study  for  both  full-CL  and 
2/3-CL  cloud  cases  for  16  January  1974,  and  comparative 
values  after  annual  climatology  of  Budyko  (1956).   Flux 
values  in  ly  min~ 1 . 
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H.   COMPARISONS  WITH  CLIMATOLOGICAL  RESULTS  OF  F1Q* 

Table  III  (b)  depicts  the  zonally-averaged  values  of 
the  Flf*  at  the  earth's  surface  obtained  in  this  study  using 
both  full-cloud  model  and  the  two-thirds  cloud  model.   A 
comparison  of  the  model -computed  F-.Q*-values  for  the  full 
and  two-thirds  cloud  cases  is  then  made  with  annual  mean 
values  of  F-0*  obtained  from  climatological  studies  (Budyko, 
1956).   The  F..  *  values  shown  in  Table  III  (b)  were  computed 
from  Equation  3-9  using  the  full  CL.. jCL-  values  or  the  re- 
duced (2/3CL..  ,  2/3CL-)  in  the  respective  cloud-models. 
Since  Budyko '  s  mean  values  of  F-in*  were  based  on  annual  cli- 
matology, his  overall  latitudinal  range  of  F..  *  values  is 
not  as  great  as  that  of  F10*  by  either  cloud-model  presented 
in  Table  III  (b)  for  16  January  1974. 

It  is  to  be  noted  that  the  locations  of  highest  values 
of  Flf*  after  Budyko  agree  with  those  computed  for  both  F-0* 
(CL)  and  F..  *(2/3  CL)  ,  namely  in  the  subtropics  (latitude 

20°N) ,  and  again  in  the  latitude  range  35-40°N.   The  "hemi- 

i 


spheric"  weighted  mean  of  F..  *  =  .096  ly   min   by  the  re- 
duced-cloud model  is  slightly  greater  than  the  cosine-weighted 
mean  of  Budyko  (.094  £y  min   ).   However,  the  weighted  mean 
of  Fin*(CL)  resulting  from  the  full-cloud  model  is  consider- 
ably smaller,  F. n*  (CL)  =  .075  ly   min   .   The  chief  differences 
in  F-0*(CL)  as  a  function  of  latitude  compared  with  that  of 
Budyko  is  accounted  for  by  the  marked  minimum  in  F..  *  cor- 
responding to  the  relatively  high  composite  CL-amounts  in 
low-latitudes  (from  20S  to  5N)  and  again  in  the  latitude  band 
45-55°N  where  in  both  areas  typical  minimum  values  of  Fin*(CL) 
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=  .055  ly   min    occur.   Such  a  minimum  of  F..  *  in  the  regions 
45-55°N  is  generally  associated  with  the  location  of  the 
polar-front  zone.   However,  in  this  same  region  the  2/3-CL 
model  gives  average  values  of  F1  *  =  .068  ly   min   .   The 
zonal-mean  values  of  CL  and  2/3-CL  applicable  to  Table  III 
(b)  by  latitude  are  exactly  as  listed  in  Table  III  (a).   It 
is  not  clear  from  the  comparisons  of  Table  III  (b)  whether 
the  higher  values  of  F..  *  deduced  from  the  2/3-CL  model 
represents  a  closer  approximation  to  reality  than  the  F..  *(CL) 
values  resulting  from  Equation  3-9  with  the  full-cloud  model 
values . 

It  is  to  be  expected  that  analogous  comparisons  made 
using  the  solar-insolation  model  (summarized  in  Table  IV) 
will  indicate  some  relative  advantage  for  using  the  full  or 
two-thirds  cloud  model. 
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IV.   SOLAR  RADIATION 

A.   PARTITION  OF  SOLAR  INSOLATION 

The  solar  constant  assumed  in  this  study  at  level  k  =  0 
(top  of  atmosphere)  was  2.00  ly   min    (Joseph,  1971).   This 
flux-value  was  further  depleted  by  4%  to  account  for  the 
attenuation  caused  by  oxygen  and  ozone  above  the  tropopause. 
This  left  the  value  S  =  1.92  ly   min   at  level  k  =  2  to  be 
used  in  this  study  as  the  effective  solar  constant. 

Equation  (4-1)  was  then  utilized  to  compute  the  effective 
solar  insolation  at  the  tropopause  (k  =  2)  as  follows: 

F(2)  =  S  [  f-  ]"2  Cos  z  (4-1) 

m 

where  S  =  effective  solar  constant  at  k  =  2 

Cos  z  =  cosine  of  the  zenith  angle  for  the  Julian  date  used 

r/r   =  ratio  of  the  actual  earth-sun  distance  to  the  mean 
earth-sun  distance  for  the  Julian  date  used  in 
this  study. 

The  Smithsonian  Meteorological  Tables  (List,  1958)  gives  the 

ratio  r/r   and  the  solar  declination  6  for  16  January,  0000GMT, 

m  '  ' 

as  these  values  are  employed.   These  values  are,  respectively: 

—  =  .983715 
rm 

6  =  -21.08  degrees  of  lat. 

6  is  used  in  evaluating  the  cosine  of  the  solar  zenith  angle, 
given  by 

Cos  z  =  Sin  <t>  Sin  6  +  Cos  <j>  Cos  6  Cos  h  (4-2) 

where  §    is  the  latitude  and  h  is  the  hour  angle  of  the  sun 
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relative  to  the  meridional  data  lines.   For  example,  Fig.  1 
makes  it  clear  that  for 

line  1         h  =  55° 

line  2         h  =  10° 

lines  3,  4     h  =  35° 
at  the  times  of  the  synoptic  charts  (0000GMT  and  1200GMT, 
respectively)  .   Sin  4>  was  computed  according  to  the  standard 
polar  stereographic  projection  formula  applicable  to  the 
FNWC  base  chart  given  by 

rF2  -  [(I-32)2  +  [(J-32)2] 

Sin  <f>  =  — (4-3) 

rE2  f  [(I-32)2  +  [(J-32)2] 

where  rp2  =  973.752.   Thus  Sin  <f>  assumes  the  following  func- 
tional form  in  terms  of  the  FNWC  grid-coordinate  I  (Fig.  1): 

Lines  1,  3,  4 

c.   .    975. 752-2  (52-1) 2  ,.     .     ,    ,, 

Sin  *  =  973.752^2(32-1)2  ^4"4  M) 

Line  2 

c.   .    973.752-(32-I)2  ,A     .     r,.. 

Sin  *  =  973.752+(52-I)2  (4"4  (b)) 

or,  conversely  I  is  given  in  terms  of  cj>  by 

Lines  1,  3,  4       I  =  32  -  22.065[1  ^°|.*  y]    (4-4  (c)) 

Line  2  I  =  32  -  31.205  [^-^sJL-]    (4-4  (d)) 

1,...,  25  for  Lines  1,  2 
I  =     9 , . . .  ,  25  for  Line  3 
63, ... ,  38  for  Line  4 . 

For  lines  1,  2,  and  3  the  gridpoint  soundings  correspond  to 

+•  \\ 
0000GMT  16  January  when  solar  noon  occurs  at  the  180 
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meridian.   Line  4  gridpoint  soundings  correspond  to  1200  GMT, 
16  January,  when  solar  noon  was  at  the  Greenwich  meridian. 

A  very  simple  partition  of  solar  insolation  was  utilized 
in  this  study  after  Joseph  [1971] .   It  consisted  of  dividing 
the  insolation  F(2)  into  two  parts  at  level  k  =  2.   One  part 
was  considered  to  include  all  wavelengths  A  >  .  9  ym  where 
absorption  by  water  vapor  and  carbon  dioxide  bands  are  the 
most  prevalent  attenuation  processes  in  clear  air.   This  part 
of  solar  insolation  was  listed  as  F(A)  energy  and  considered 
subject  to  water-vapor  absorption  but  not  to  Rayleigh  scat- 
tering.  For  those  wavelengths  X    <  .9  ym,  absorption  of  the 
solar  insolation  energy  by  water  vapor  was  considered  neg- 
ligible.  This  part  of  the  solar  insolation  band  was  listed 
as  F(S)  and  was  subject  only  to  Rayleigh  scattering  attenua- 
tion in  clear  air.   The  two  partitions  are  formulated  as 
follows : 

F(A)  =  .349  F(2)  (4-5  (a)) 

F(S)  =  .651  F(2) .  (4-5  (b)) 

In  this  study,  the  introduction  of  two  cloud  decks  produced 
cloud-reflectivity  effects  upon  both  the  F (A)  and  F(S)  solar 
energy  insolations.   However,  in  the  clear  areas  around  any 
gridpoint  the  absorption-attenuation  only  applies  to  the 
F(A)  insolation  while  the  Rayleigh  scattering-attenuation 
applies  to  the  F(S)  insolation  only. 

B.   DISPOSITION  OF  F(S)  INSOLATION 

In  the  treatment  of  the  F(S)  insolation,  Joseph  [1971] 
determined  that  Rayleigh  scattering  reflectance  in  clear  skies 
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(after  Coulson,  1959)  could  be  effectively  approximated  by 
least  squares  in  the  following  form 

a(R)  =  .085  +  .25074  [log  (  J-  Sec  z  )]        (4-6) 

o 

where  P   =  1013.25  mb .   In  Equation  4-6,  tt/P   =  1  in  view  of 
o  o 

the  fact  that  mean  sea  level  pressure  tt  is  close  to  1000  mb . 
Also 

Sec  z  =  (Cos  z) 

with  Cos  z  given  by  Equation  4-2. 

The  surface  albedo  cx(G)  is  another  reflective  parameter 
utilized  in  this  study.   Over  oceanic  areas  the  following 
formula  for  a(G)  after  Gates  et   al .  [1971]  ,  was  utilized: 

oc(G)  =  max  {.06,  .06  +  .54  (.7  -  Cos  z)  }  .      (4-7) 

1 .  Clear  Sky  Case 

In  the  clear  sky  (0,  0)  case  the  F(S)  insolation  was 
subjected  only  to  Rayleigh  scattering  reflectance  a(R)  and 
to  surface  reflectance  a(G) .   Considering  the  likelihood  of 
a  succession  of  multiple  reflections  between  earth  and 
atmosphere,  the  F(S)  insolation  actually  penetrating  the 
earth's  surface  after  scattering  is  given  by 

IS10(0,0)  =  F(S)  [l-a(R)]  [  1  +  a  (R)  ct  (G)  +  .  .  .  (a(R)a(G))n 

+. ..]*(l-a(G))  (4-8  (a)) 

that  is,  by 

IS10(0,0)  =  F(S)  [l-a(R)]  [  1-a  (G)]/ [1-a  (R)  a(G)   (4-8  (b)) 

2 .  Cloudy-Sky  Cases 

In  the  three  cases  in  which  clouds  were  present,  F(S) 
insolation  absorbed  by  the  ground  at  each  gridpoint  was 
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computed  using  a  variation  of  the  following  equation  (after 
Arakawa,  1972) : 

IS10(1,1)  =  F(S)(l-R(l))(l-R(2))(l-a(G)) 

*{l-[(R(l)R(2)+R(2)a(G)+R(l)a(G) 

+  2R(l)R(2)a(G)]}"\  (4-9) 

As  indicated  by  the  notation  (1,  1),  denoting  CL(1)  =  CL(2)  = 
1.0.   Equation  4-9  is  the  formula  used  in  calculating  F(S) 
insolation  absorbed  by  the  ground  in  the  case  where  overcast 
clouds  are  present  at  both  levels  of  Fig.  2.   Also  in  Equa- 
tion 4-9,  constant  cloud-reflectance  values  were  chosen, 
namely  R(l)  =  .54  for  the  mid-level  clouds  between  k  =  4 
and  6,  and  R(2)  =  .66  for  the  low-level  clouds  between  k  = 
8  and  9.   Both  cloud-reflectance  values  are  as  suggested  by 
C.  D.  Rodgers  [1967] . 

For  the  other  cloud  cases,  the  following  changes 
were  applied  to  Equation  4-9.   In  the  (1,  0)  case  (CL(1)  = 
1.0,  CL(2)  =  0.0),  consider  that  R(2)  =  0  in  (4-9),  from 
which  it  follows  that 

IS10(1,  0)  =  F(S)(l-R(l))(l-a(G))/[l-R(l)a(G))]. 

(4-10) 

In  the  case  (0,  1),  R(l)  =  0.0  in  (4-9)  so  that 

IS10(0,  1)  =  F(S)(l-R(2))(l-a(G))/[l-R(2)a(G)].(4_ii) 
3 .   Composite  F(S)  Insolation 

Equations  (4-8),  (4-9),  (4-10),  and  (4-11)  were 
utilized  in  the  computation  of  the  cloud-weighted  F(S)  in- 
solation penetrating  the  earth's  surface  considering  the 
areal-weights  of  the  cloud  combinations  denoted  by  (0,  0), 
(1,  1),  (1,  0)  and  (0,  1)  about  a  gridpoint.   The  resultant 
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F(S)  insolation  penetrating  the  earth's  surface  denoted  by 
IS10  is  therefore  expressible  as 


(4-12) 


IS10(CL(1),  CL(2))  =  IS10(0,  0)  W(0,  0) 

+  IS10(1,  1)  W(l,  1) 
+  IS10(1,  0)  W(l,  0) 
+  IS10(0,  1)  W(0,  1) . 

Here  the  weighting  factors  W(0,  0),  W(l,  1),  W(l,  0)  and 
W(0,  1)  are  as  computed  in  Equations  3-1,  3-2,  3-3,  and  3-4 
respectively.   Note  finally  that  the  part  of  F(S)  insolation 
reflected  to  space  is  found  by  subtracting  IS10(CL(1),  CL(2)) 
from  F (S) . 

C.   DISPOSITION  OF  F (A)  INSOLATION 

The  fractional  portion  of  the  solar  insolation  subject 
to  absorption  by  atmospheric  water-vapor  and  carbon  dioxide 
are  covered  in  the  following  subsections. 

1.   The  Clear-Sky  Case  (0,  0) 

The  Manabe-Moller  absorptivity  function  provided  the 
necessary  absorptivity  values  for  the  key  layers  in  this  case 
The  form  of  this  absorptivity  function  is 

a(2,  k)  =  .271[U(2,  k)Sec  z]'303.  (4-13) 

Here  a_  is  the  absorptivity  applied  to  the  pressure-scaled 
water  vapor  mass  between  levels  2  and  k  (Fig.  2)  along  the 
zenith  slant-path  angle  z.   The  resultant  absorbed  insola- 
tional  energy  in  the  particular  layer  (2,  k)  is  then  given 
by  the  Manabe-Moller  relation 

A(2,  k)  =  0.271F(A) [U(2,  k)  Sec  z]*303.       (4-14) 
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The  two  layers  of  interest  in  which  absorption  was  computed 
were  (2,  6)  and  (2,  10).   The  absorbed  insolation  in  the 
layer  (6,10)  was  then  computed  by 

A(6,  10)  =  A(2,  10)  -  A(2,  6).  (4-15) 

It  is  noted  here  that  the  water-vapor  mass  above  level  2  was 
assumed  negligible  in  this  study. 

By  subtracting  A(2,  10)  from  F  (A) ,  the  direct  transmis- 
sion of  F (A)  insolation  impinging  at  the  earth's  surface  was 
determined.   The  transmission  of  F(A)  insolation  is  then  fur- 
ther reduced  by  the  transmissivity  (l-a(G)),  after  surface- 
reflectance  which  leads  to  the  earth-absorbed  result 

IA10(0,  0)  =  F(A){l-.271[U(2,10)Sec  z] * 303} (1-a (G) ) . 

CM-16  (a)) 

The  transmitted  energy  impinging  upon  the  earth  just  prior  to 
absorption  is 

TRANA(0,  0)  =  IA10(0,  0)  /  [l-a(G)].  (4-16  (b)) 

2 .   Cloudy  Cases 

In  order  to  compute  meaningful  dispositions  of  F(A) 
insolation  in  cloudy-sky  cases,  cloud  reflectivities  and 
cloud  absorptivities  after  C.  D.  Rodgers  [1967]  were  utilized. 
The  reflectance -values  used  here  are  different  from  those 
suggested  by  Rodgers  for  the  F(S)  wavelengths.   The  cloud 
reflectivities  used  here  are  RA(1)  =  .46  and  RA(2)  =  .50. 
In  this  case  there  are  also  cloud-absorptivities  to  be  con- 
sidered.  These  were  taken  as  A(l)  =  .20  and  A(2)  =  .30, 
respectively.   In  the  following  discussions  the  cloud  con- 
ditions are  considered  totally  overcast  in  each  of  layers 
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indicated  when  the  notations  ((1,  1),  (0,  1),  and  (1,  0)) 
are  utilized. 

Schematic  representations  of  the  computations  per- 
formed in  the  various  cases  ((1,  1),  (0,  1),  (1,  0))  are 
displayed  in  Figures  7,  8,  and  9.   These  figures  suggest  the 
parameters  required  in  computing  the  disposition  of  incoming 
F(A) -insolation  from  level  k  =  2  to  the  earth's  surface 
(k  =  10). 

Figure  7  depicts  the  disposition  of  F(A)  insolation 
in  the  overcast  case  with  both  high  and  low  cloud  layers 
(1,  1).   Figure  8  shows  the  case  (1,  0)  with'  an  upper  over- 
cast only.   Figure  9  displays  the  case  (0,  1)  for  a  low 
overcast.   The  equations  (after  Warner,  1974)  which  relate 
to  the  parameters  in  these  figures  are  listed  in  Appendix 
B-l,  B-2,  B-3,  respectively. 

In  Figs.  7,  8,  and  9,  the  insolations,  A24,  A46, 
A68,  A89,  and  A910,  etc.,  represent  the  insolation  absorbed 
in  the  layers  concerned.   Symbols  F2,  F4,  F6,  F8,  and  F9, 
etci,  depict  the  streams  of  insolation  passing  through  the 
indicated  level.   A  vertical  arrow  implies  the  direction  of 
insolation  passage,  i.e.,  4-  denotes  downward  insolation,   I 
upward-reflected  insolation,  and  4-4-  downward-reflected  inso- 
lation.  The  absorption  quantity  A(6,  8)4-4-,  for  example,  in- 
dicates absorbed  energy  remaining  in  (6,  8)  from  a  downward 
reflected  beam. 

Since  multiple  reflections  occur  between  the  earth's 
surface  and  a  cloud  base,  or  between  cloud-layers  in  this 
model,  the  insolational  amount  reaching  the  lowermost  reflecting 
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Figure  7    Schematic  representation  of  F(A)  insolation 
disposition  in  the  case  of  two  overcast  layers. 
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Figure  8    Schematic  representation  at  F(A)  insolation  dis- 
position with  an  upper  overcast  layer  only. 
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Figure  9    Schematic  representation  of  F(A)  insolation 
disposition  with  a  lower  overcast  layer  only. 
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surface  is  affected  by  the  downward  reflected  beams.   Addi- 
tional insolational  amounts  remaining  after  more  than  two 
reflections  are  negligible  and  are  not  included  in  the  models 
(Appendix  B) .   Also  note  that  insolation  reflected  upward 
from  a  lower  interface,  either  cloud-top  or  ground,  to  the 
base  of  a  higher  cloud  deck  was  not  subjected  to  absorption 
by  or  transmission  through  the  cloud.   This  simplification 
resulted  in  slightly  reduced  F (A)  insolation-reflectance  to 
space . 

A  final  point  to  note  is  that  there  are  three  con- 
tributions to  the  absorption  of  solar  insolation  within  a 
layer  between  any  two  reflecting  surfaces.   In  the  Figs.  7, 
8,  9,  these  contributions  are  identified  by  arrows  which  in- 
dicate the  portion  of  path  being  crossed. 

For  each  of  the  overcast-cloud  combinations,  the 
transmitted  F(A)  insolation  arriving  at  the  earth's  surface 
may  be  defined  using  the  following  notation  (TRANA)  as 

TRANA(1,  1)  =  F10I  +  F104-4-  (4-17) 

TRANA(1,  0)  =  F104-  +  F1(H+  (4-18) 

TRANA(0,  1)  =  F10I  +  F10-K  (4-19) 

where  the  right  side  parameters  of  (4-17),  (4-18),  and  (4-19) 
are  derived  in  Appendices  (B-l),  (B-2),  and  (B-3) . 

The  F(A)  insolation  absorbed  by  the  earth  in  each 
case  is  derived  by 

IA10(1,  1)  =  F1(H  (l-a(G))  +  F1(H+  (4-20) 

IA10(1,  0)  =  F10I  (l-a(G))  +  F104-4-  (4-21) 

IA10(0,  1)  =  F10!  (l-a(G))  +  F10+1.  (4-22) 
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In  the  last  three  equations,  the  quantity  F104-+  is  small 
enough  in  each  case,  that  no  further  reflections  from  the 
earth  were  considered. 

3 .  Composite  F (A)  Layer-Absorptions  and  Surface- 
Absorption  Insolation 

As  has  been  previously  discussed,  the  standard  grid- 
area  weighting  scheme  of  this  study  was  applied  to  obtain 
composite  values  of  the  absorbed  F(A)  insolation  in  key 
layers  and  also  within  the  earth's  surface.   The  weighting 
factors  applied  to  the  corresponding  overcast-combination 
absorption  quantities  provided  the  following  composite 
results : 

A26(CL(1),  CL(2))  =  A26,0   0^W(0,  0)  +  kit*  ^      Q)W(1,  0) 

+  A26(0^  1}W(0,  1)        (4-23) 
+  A26(1>  1}W(1,  1) 

A610(CL(1),  CL(2))  =  A610(0   0)W(0,  0) 

+  A610(lj  0)W(1,  0)       (4-24) 

+  kblO^      X)W(1,  1)+  A6l0(Q   1)W(0J  1) 

IA10(CL(1),  CL(2))  =  IA10(0,  0)W(0,  0) 

+  IA10(1,  0)W(1,  0) 

(4-25) 
+  IA10(0,  1)W(0,  1) 

+  IA10(1,  1)W(1,  1) . 

The  weighting  factors  W(0,  0),...,W(1,  1)  were  first  listed 

in  Equations  3-1,  3-2,  3-3,  3-4. 

4.  Absorptivity  (ABA)  by  Layers 

Here  the  (fractional)  absorptivity  as  well  as  the 
actual  insolation  values  absorbed  in  the  layers  are  considered. 
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In  the  computation  of  absorptivity,  the  total  undepleted  in- 
solation at  the  top  (k  =  0)  is  used.  The  following  equation 
was  utilized  in  this  calculation: 

FADJ  =  2.00(r/r  )~2  Cos  z.  (4-26) 

ABA,  or  the  absorptivity  of  the  troposphere  was  computed  from 
the  ratio  of  the  insolation  absorbed  in  the  troposphere  to 
the  insolation  incident  at  the  top  of  the  atmosphere  rather 
than  at  k  =  2 : 

ABA  =  A(2>  6)FADj(6>  1Q)  •        •  (4-27) 

D.   ALBEDO  (ALB)  OF  THE  EARTH-TROPOSPHERE  SYSTEM 

In  considering  the  planetary  albedo,  the  absorptions  of 
the  earth-troposphere  system  in  both  the  F(A)  and  F(S)  in- 
solational  regions  must  be  recalled  by  the  program.   Thus 
we  define  the  reflected  insolational  energy  in  F(A)  as  REFA, 
and  it  is  computed  by 

REFA  =  F(A)  -  A26  -  A610  -  IA10.  (4-28) 

Likewise  REFS,  the  reflected  part  of  F(S)  is  defined  by 

REFS  =  F(S)  -  IS10.  (4-29) 

It  is  understood  that  REFA,  REFS  are  computed  for  appropriate 
CL,,  CL2  at  each  gridpoint  sounding  considered. 
Summing  REFA  and  REFS  gives  REF: 

REF  =  REFS  +  REFA.  (4-30) 

Finally  the  planetary  albedo  is  related  to  FADJ  through 

ALB  ■  mn  •  C4"3" 
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E.   COMPOSITE  ABSORPTIVITY  (ABG)  BY  THE  EARTH-SURFACE; 
COMPOSITE  ATMOSPHERIC  TRANSMISSIVITY  (ATRAN) 

1.  Absorptivity  (ABG)  of  Earth 

By  summing  the  weighted  values  of  F(S)  and  F(A)  por- 
tions of  the  incoming  insolation  entering  the  earth,  the 
total  insolation  absorbed  at  the  earth's  surface  was  com- 
puted.  This  quantity  when  divided  by  the  extraterrestrial 
insolation  gave  the  fractional  absorptivity  (ABG)  of  the 
earth's  surface.   The  equation  for  ABG  was 

AD~    IA10  +  IS10  fA    rn, 

ABG  =  FADJ (4"32) 

where  IA10,  IS10,  and  FADJ  were  defined  previously  by  Equa- 
tions 4-25,  4-12,  and  4-26,  respectively. 

2 .  Transmissivity  (ATRAN)  of  the  Troposphere 

Also  computed  was  the  total  insolational  energy 
TRAM,  incident  at  the  earth's  surface  just  before  absorption 
by  the  surface.   This  calculation  is  given  by 

TRAN  =  TRANA  +  [IS10  /  (l-a(G)].  (4-33) 

Here  a(G)  is  as  previously  defined  in  Equation  4-7.   TRANA 
is  as  defined  by  the  weighted  value  of  TRANA(0,  0),  TRANA (1,  1) , 
TRANA (1,  0)  and  TRANA (0,  1)  given  by  Equations  4-16,  4-17, 
4-18,  and  4-19.   Also  note  that  the  four  terms  of  IS10  of 
(4-8),  (4-9),  (4-10),  and  (4-11)  have  the  common  factor 
(l-a(G))  in  the  numerator  and  therefore  each  F(S)  insolation 
component  at  the  earth  just  before  absorption  need  only  be 
divided  by  (l-a(G)).   TRAN  may  thus  be  viewed  as  the  total 
insolational  energy  incident  at  a  pyrheliometer  located  at 
earth.   The  (fractional)  transmissivity  of  the  troposphere 
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(ATRAN)  is  then  computed  from 

ATRAN  E  TRAN  /  FAD J .  (4-34) 

Note  finally  that  the  major  dispositions  of  the  total  in- 
solation at  the  indicated  map  times  have  been  identified  by 
the  fractional  values,  ALB,  ABA,  ABG,  and  ATRAN,  each  of 
which  is  a  fractional  value  representing  the  reflectivity 
(albedo),  absorptivity,  or  atmospheric  transmissivity  as 
the  case  may  be. 

3 .   Computational  Check 

The  computational  scheme  utilized  in  this  model  was 
checked  by  summing  the  fractional  values  ALB,  ABA,  and  ABG 
relative  to  the  troposphere  at  each  gridpoint.   The  value  in 
each  case  must  total  .96,  since  as  previously  noted  the  at- 
tenuation of  solar  insolation  was  taken  as  four  percent  as 
it  streamed  through  the  stratosphere. 

F.   STATISTICAL  ANALYSIS 

In  order  to  substantiate  some  of  the  computations  per- 
formed in  this  section,  several  of  the  most  important  items 
computed  were  statistically  analyzed  using  linear  regression 
computer  programs  from  the  BMD  set  of  statistical  programs 
(Dixon,  1973). 

1 .   Clear  Sky  Cases 

Using  ALB  (0,  0,  z) ,  which  is  the  clear  sky  case  of 
albedo,  as  the  predictand,  and  log,pSec  z  and  its  square  as 
the  predictor,  the  following  best-fit  equation  resulted 

ALB(0,  0,  z)  =  .05705  +  .54529  loginSec  z 

iU         (4-35) 
-  .18885(log  Sec  z)2. 
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R   =  .9918 
m 

where  R  represents  the  multiple  correlation  coefficient. 

This  result  was  as  expected  over  the  ocean  where 

a(G)  and  a(R)  involved  logarithmic  dependence  on  Sec  z. 

Average  values  of  ALB  and  logir|Sec  z  were 

ALB(0,0)  =  .20350  log^Sec  z  =  .34049. 

Other  clear  sky  regression  tests  made  use  of  the  water-vapor 
mass  path  length  (M)  defined  by 

M  =  (U  Sec  z)h    log10(U  Sec  z)^.  (4-36) 

This  parameterization  of  water  vapor  mass  is  similar  to  that 
developed  by  Hanson  [1971],   who  used  a  similar  M  as  a  pre- 
dictor in  his  empirical  formulations  of  ABA  for  both  clear 
and  partly  cloudy  sky  cases.   In  this  study  the  two  regres- 
sions attempted  using  M  as  a  predictor  were  with  ABA(0,  0,  M) 
and  ATRAN(0,  0,  M)  as  predictands.   The  best-fit  equations 
resulting  were 

ABA(0,  0,  M)  =  .11478  +  .04012  M  (4-37) 

R   =  .9675 
m 

ATRAN(0,  0,  M)  =  .77794  -  .06337  M  (4-38) 

R      =  .8575. 
m 

The  means  of  ABA(0,  0,  M)  and  ATRAN(0,  0,  M)  for  the  original 

93  soundings  were 


ABA(0,  0,  M)  =  .15154 


ATRAN(0,  0,  M)  =  .71987 

M  =  1.83259  (gm  cm"2) 
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2 .   Relationships  Between  Albedo,  Atmospheric  Absorptivity, 
Atmospheric  Transmissiyi ty ,  and  Ground  Insolation  in 
the  Cloudy  and  Clear  Sky  Cases 

In  addition  to  the  clear  sky  value  of  ALB,  it  was 

also  possible  to  compute  the  composite  cloudy  sky  case.   A 

regression  was  formed  showing  the  relationship  between  the 

ratio  ALB(CL(1) ,CL(2) ) /ALB (0 ,  0)  as  the  predictand  and  the 

total  opaque  cloud  cover  (CL)  as  the  predictor.   CL  in  this 

model  is  specified  by  Equation  4-39  below: 

CL  =  CL(1)  +  CL(2)  -  CL(1)CL(2).  (4-39) 

CL ,  it  was  felt,  gave  a  good  approximation  of  the  effective 
cloud-cover  by  the  two  layers  of  cloud  amounts  CL.  and  CL? 
by  Equation  2-14  used  in  this  model.   As  explained  earlier 
in  this  study,  a  one-third  reduction  of  both  CL(1)  and  CL(2) 
after  initial  determination  by  Equation  2-14  was  also  tried 
in  all  computations  involving  cloud  amounts.   This  test 
represented  an  initial  attempt  at  tuning  the  cloud  model  for 
radiative  calculations  as  will  later  be  substantiated,  and 
gave  results  in  closer  agreement  with  the  latest  literature 
on  satellite  reflectances.   Thus  all  subsequent  regressions 
involving  CL  will  also  include  those  best-fit  equations  for- 
mulated using  a  CL  based  also  upon  a  one-third  reduction  of 
CL(1)  and  CL(2).   A  small  letter  subscript  'a*  in  the  equa- 
tion number  will  indicate  those  computations  using  full-cloud 
amounts  while  a  small  letter  subscript  Tb'  will  indicate  those 
utilizing  a  one-third  reduction  in  the  cloud  model  amounts. 

The  first  regression  equations  tested  were  those  for 
ALB(CL)  and  ABA(CL)  for  the  respective  cases  of  "full"  and 
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two-thirds  cloud  cover.   The  results  are  given  in  (4-40) , 
4-41) .   The  symbol  R   once  again  signifies  the  multiple 
correlation  of  the  statistical  regressions. 

ALB  (CL(1)CL(2)  =  ALB(0,0)[1  +  3.59683CL 
a 

-  1.7006CL2]  (4-40  (a)) 

R  =  .8755 
m 

ALBb(CL(l) ,CL(2)  =  ALB(0,  0)[1  +  4.13667CL 

-  3.18972CL2]  (4-40  (b)) 

R      =  .8876. 
m 

The  sample  means  of  the  values  in  Equations  4-40a  and  4-40b 

are 

Equation  4-40a  Equation  4-40b 

ALB(CL(1)  ,CL(2))  =  .4298        ALB(CLO-)  ,CL(2))  =  .3640 


ALB(0,  0)         =  .2035        ALB(0,  0)         =  .2035 
CL  =  .6951        CL  =  .5176 

The  correlation  coefficient  in  Equation  4-40  appears  to  show 
strong  dependence  of  ALB  upon  CL.   Also  note  that  ALB(CL(1), 
CL(2))  values  computed  here  are  slightly  higher  than  reported 
by  Raschke  et    al. 3     [1973],  (ALB  =  .30)  for  essentially  the 
same  set  of  gridpoints  during  the  NIMBUS  III  period  21  January 
through  3  February  1970.   However,  the  one-third  reduction-CL 
appears  to  give  much  closer  albedo  results,  suggesting  that 
the  full  cloud  amounts  of  Smagorinsky  by  2-14  are  too  high. 

A  second  regression  computed  was  between  the  ratio  of 
tropospheric  absorptivity  of  the  cloudy  air  to  that  for  the 
clear  sky  case,  versus  CL.   The  following  best-fit  equation 

resulted: 
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ABA  (CL(1),CL(2))  =  ABA(0,  0)[1  +  .5043CL 

a  ■  (4-41  (a)) 

-  .1326CL2] 

R  =  .8989 
m 

ABA,(CL(1),CL(2))  =  ABA(0,  0)[1  +  .5522CL 

(4-41  (b)) 

-  .2451CL2] 

R   =  .8873 
m 

Mean  statistics  for  this  regression  case  are 

Equation  4-41(a)  Equation  4-41(b) 

ABA  (CL(1) ,CL(2))  =  .1916       ABA,(CL(1) ,CL(2j  =  .1810 
a  d 

ABA(0,  0)  =  .1515       ABA(0,  0)         =  .1515 


CL  =  .6951       CL,  =  .5176 

a  b 

In  (4-41),  it  is  seen  that  the  model  specifies  an  increase  in 

atmospheric  (solar)  absorptivity  with  increasing  the  cloud 

cover.   This  result  is  in  agreement  with  studies  of  Plante 

[1972],  Warner  [1974],  and  with  the  earlier  study  of  London 

[1957]  . 

An  analogous  statistical  regression  was  then  developed 

for  the  cloudy-sky  transmissivity  relative  to  the  clear-sky 

transmissivity .   The  resulting  regressions  were 

ATRAN  (CL(1) ,CL(2))  =  ATRAN (0 , 0) [ 1- . 6020CL 

E  (4-42  (a)) 

-  .0285CL2] 

R   =  .9895 
m 

ATRAN,  (CL(1),CL(2))  =ATRAN(0,0)  [1-.6833CL 

(4-42  (b)) 
+  .1246CL2] 

R   =-  .9862. 
m 
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Equation  4-42a  Equation  4-42b 

ATRAN  CL(1) ,CL(2)  =  .4440   ATRAN,  CL(1) ,CL(2)  =  .4957 
a  d 

ATRAN(0,  0)       =  .7199   ATRAN(0,  0)       =  .7199 

If  Equation  4-42b  is  terminated  at  the  first  power  of  CL,  the 
result  reduces  to 

ATRAN (CL(1) ,CL(2))  =ATRAN(0,  0) [1- . 5961CL]   (4-43) 

R   =  .9855. 
m 

The  result  (4-43)  is  in  reasonable  agreement  with  the  Savino- 

Angstrom  result  which  is 

ATRAN(CL(1) ,CL(2))  =  ATRAN (0,  0) [1- (1-k) CL]   (4-44) 

Here  k  is  a  slowly  increasing  function  of  latitude,  ranging 
from  0.35  to  0.50,  according  to  Budyko  [1956]. 
3 .   Some  Conclusions 

The  statistical  analyses  performed  seem  to  show  agree- 
ment with  observational  results  of  other  investigators.   Thus 
the  radiation  model,  with  the  solar  seasonal  effects  imposed 
appears  to  be  realistic.   Further,  the  one-third  reduction  in 
amounts  given  by  the  cloud  parameterization  model  appear  to 
give  even  better  agreement  with  the  latest  satellite  findings 

G.   ALEEDO  COMPARISONS  WITH  PUBLISHED  RESULTS 

The  tropospheric  albedo  computations  of  the  solar-inso- 
lation model  of  this  section  have  been  presented  for  both 
the  full-cloud  and  two-thirds  CL  cases,  respectively.   These 
albedos  have  been  interpolated  to  whole  multiples  of  5° 
latitude  between  20°S  and  65°N.   These  computations  were  made 
for  each  of  the  four  oceanic  meridians  and  the  resulting 
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zonally-averaged  albedos  are  presented  as  a  function  of 
latitude  in  Table  IV  for  both  the  full-CL  and  two-thirds  CL 
cases,  where  comparison  is  also  made  with  satellite  clima- 
tology of  Raschke  et   al .  [1973] ,  pertaining  to  the  NIMBUS 
III  averaged  albedo  measurements  over  the  period  21  January  - 
3  February  1970.   For  sake  of  consistency  the  Raschke  albedo 
analyses  have  been  derived  by  interpolation  to  latitude  using 
the  same  four  meridians  before  applying  the  zonal-averaging 
process . 

The  model  results  under  the  headings  "FULL-CL"  and  "2/3- 
CL"  correspond  to  the  composite  cloud  amounts  previously 
listed  in  Table  11(a)  according  to  the  equation 

CL  =  CL1  +  CL2  -  CL1  *  CL2  (3-17) 

with  either  (CL-,  CL2)  by  Equation  2-14  or  (CL'   CL')  by 
Equation  2-15. 

The  most  unusual  results  to  be  observed  in  Table  IV  are 
the  small  albedo-values  (~  .2)  between  20°S  to  20°N  after 
Raschke,  as  contrasted  with  the  corresponding  values  (in  the 
range  . 3  to  .4)  by  2/3-CL  cloud  model  ,  and  more  markedly  with 
those  of  the  FULL-CL  model  (values  in  the  range  .35  to  .5). 
The  "hemispheric"  cosine  weighted  mean  albedos  are  as  follows: 


FULL-CL,  ALB  -  .46 


2/3-CL,  ALB   =  .38 


RASCHKE-ALB   =  .30 

At  every  latitude  (<J)  <  50°)  under  comparison  the  2/3-CL 
albedo  computation  is  closer  to  Raschke 's  observations  than 
is  the  full-CL  value.   At  latitudes  <J>  >  55°,  there  is  evidence 
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ALBEDO 

ALBEDO 

ALBEDO 

FULL-CL 

(2/3-CL) 

Lat 

(FULL-CL) 

(2/3-CL) 

(Raschke) 

AMOUNT 

AMOUNT 

20S 

.5462 

.4458 

.2195 

.932 

.758 

15 

.5380 

.4336 

.  2050 

.927 

.733 

10 

.4901 

.3859 

.1875 

.867 

.634 

5 

.4334 

.3432 

.1967 

.746 

.536 

0 

.4447 

.3478 

.2325 

.790 

.556 

5 

.4290 

.3362 

.2225 

.766 

.527 

10 

.3532 

.2859 

.2050 

.557 

.376 

15 

.3517 

.2895 

.2075 

.532 

.342 

20 

.3390 

.2868 

.2450 

.443 

.279 

25 

.3664 

.3151 

.2925 

.463 

.320 

30 

.4217 

.3668 

.2725 

.559 

.411 

35 

.4145 

.3703 

.3375 

.516 

.382 

40 

.4595 

.4275 

.3925 

.628 

.529 

45 

.5558 

.5151 

.4250 

.833 

.740 

50 

.5682 

.5186 

.5325 

.911 

.759 

55 

.5717 

.5218 

' .5929 

.931 

.749 

60 

.5180 

.4848 

.5600 

.764 

.528 

65N 

.4294 

.4153 

.7000 

0.0 

0.0 

Wt. 
Avg 

.4622 

.3820 

.3010 

.688 

.514 

Table  IV.   Comparison  of  planetary  albedo  as  found  by  this 
study  for  both  full  and  two-thirds  cloud  models  for  16 
January  1974  and  by  Raschke  et  al.  (1973)  based  upon  NIMBUS 
III  measurements.   Also  included  are  the  globally  weighted 
mean  values  at  the  base  of  the  columns  and  the  composite 
cloud  amounts  for  both  cloud  cases. 
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of  added  surface-ice  reflections  contributing  to  the  Raschke 
values  of  planetary  albedo.   This  is  especially  true  at  the 
latitude  65°N  where  Raschke's  results  indicate  ALB(65°N)  = 
.7,  whereas  the  model -computations  with  essentially  zero 
cloud  cover  is  close  to  .42.   Part  of  the  discrepancy  with 
the  model  at  high  latitudes  lies  in  the  non-inclusion  of  a 
surface-ice  reflectance  parameter. 

To  summarize,  apart  from  the  region  covered  by  sea-ice, 
there  is  a  marked  tendency  for  the  2/3-CL  cloud  model  to  give 
reduced  albedos  as  suggested  by  Raschke  et   al.     [1973] ,  Von  der 
Haar  and  Hanson  [1969],  and  others.   This  result  of  the  2/3-CL 
model  is  especially  more  clearly  defined  in  tropical  latitudes 
where  Von  der  Haar  and  Hanson  [1969]  have  discussed  the  re- 
ality of  both  reduced  cloud-covers  and  the  resultant  reduced 
global-albedo.   These  latter  effects  than  must  give  rise  to 
larger  net  radiation  into  the  oceanic  surface,  which  they  also 
prescribe . 

In  Section  VI,  heat-budget  computations  along  each  of 
the  four  meridians  are  presented.   These  computations  make 
use  only  of  the  2/3-CL  cloud  model  for  purposes  of  radiative 
calculations.   This  selection  (2/3-CL)  is  supported  by  the 
results  of  Section  VII.   However,  it  is  felt  that  the  grid- 
point  method  used  here  for  2/3-CL  is  not  the  final  answer  in 
the  tropics.   In  the  tropics  the  peculiar  reflection-charac- 
teristics of  cumulus  convective  elements  must  be  more  real- 
istically accounted  for  than  by  the  large-scale  layering 
considered  in  the  solar  insolation  model  of  this  section. 
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V.   SENSIBLE  AND  LATENT  HEAT  TRANSPORT  AT  THE  SEA-AIR  INTERFACE 

A.  GENERAL 

In  order  to  compute  heat  balances  involving  the  atmosphere 
and  the  ocean-surface,  it  was  necessary  to  compute  vertical 
fluxes  of  sensible  and  of  latent  heat  across  the  air-sea  in- 
terface.  The  model  adopted  to  describe  the  turbulent-fluxes 
in  the  planetary  boundary  layer  was  essentially  that  used  in 
the  FNWC  primitive  equation  model  as  delineated  by  Kesel  and 
Winninghoff  [1972]  and  Kaitala  [1974].   The  model  as  used  in 
this  study  was  adapted  in  such  a  way  that  it  gave  similar 
quantitative  bulk-transfer  results  across  the  air-sea  inter- 
face as  the  model  used  by  Kaitala  [1974].   This  was  a  major 
deviation  from  Warner's  study  [1974],  where  he  required  a 
positive  sea-to-air  sensible  heat  transport  based  primarily 
upon  a  semi -empirical  Bowen  ratio. 

B.  SENSIBLE  HEAT  TRANSPORT 

As  previously  mentioned  a  slightly  different  modelling  as 
compared  to  that  of  Kesel  and  Winninghoff  [1972]  was  used  in 
this  study  for  the  vertical  distribution  of  sensible-heat  con- 
vergence in  the  layer  800  to  1000  mb .   In  the  presently  opera- 
tional heat-package  of  FNWC,  as  described  by  Kaitala  [1974] 
and  by  Kesel  and  Winninghoff  [1972],  the  planetary  boundary 
layer  sensible  heat  transport  is  given  by 

K*  eq  "  6x 

Hr(900)  =  p  C  (  -£ ){  Y   -  /  .  /   >  . 

P  1  +  a*  (  69  -  6g  9     10 

z9   z10  (5-1) 


8H 


In  (5-1),  the  parameters  are  as  described  below: 

C   =  .239  cal  gm"1 (deg.  K) " l 

0   =  T 
x    x 

Y   =  (the  critical  countergradient  lapse  rate)  = 

3  x  10"5(deg.  K)cm_1 
K*  =  (the  neutral  value  of  the  eddy  coefficient)  = 

2  x  105cm2sec_1 
a*  =  (a  turbulent  transfer  coefficient)  = 

5  x  10ltcm(deg.  K)  ~  * 
P10  =  (surface  density)  =  p1Q (RdT1Q) " * . 

The  other  parameters  in  (5-1)  are  0g,  which  is  the  potential 
temperature  at  900  mb ,  9   which  is  the  potential  temperature 
at  the  top  of  the  constant  flux  layer,  and  6   which  is  the 
potential  temperature  at  the  surface. 

Equation  5-1  was  given  by  Kaitala  [1974]  as  the  planetary 
boundary  layer  equivalent  of  the  bulk  formula  for  H  ,  given 
in  terms  of  surface-layer  parameters  as 

Hr  =  pinC  CnV  (T  -  T  )  (5-2) 

r    K10  p  D  s v  g   xJ  v    * 

where  Cn  is  the  drag  coefficient,  V   =  surface  wind,  and 

(T   -  T  )  is  proportional  to  the  temperature  lapse  for  the 
g    x 

thin  layer  between  the  sea  surface  and  the  top  of  the  con- 
stant-flux layer.   If  H„  of  (5-1,  5-2)  are  always  equal,  then 
it  follows  that  T   may  be  solved  from  the  identity,  giving 


-    D9 
CVT   +Kt^--KY 
t  =  P  s  g    A  Az    ^  'c 


(5-3  (a)) 

—  +  r   V 
Az    L,Dvs 
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where 

k  =  K*  6q  -  e,  t5"3  ^ 

1  +  a*  (  — s_  ) 

z9  "  z10 

Rd  9"10     10 
Az  =  zg  -  z1Q  =  —  rjr    (Jin  -g  )  (5-3  (c)) 

The  requirement  of  identical  values  of  Hf  at  level  x  and  at 
level  k  =  9,  with  H„  decreasing  to  zero  at  k  =  8,  was  inter- 
preted by  Kesel  and  Winninghoff  [1972]  as  giving  rise  to 
sensible-heat  convergence  Hr ,  layer-centered  at  k  =  9 . 
Actually  Fig.  10a  suggests  that  the  vertical  convergence  of 
sensible  heat  flux  of  amount  Hr  occurs  in  the  layer  (8,  9) 
when  the  assumption  of  constant  heat  flux  in  the  layer  (x,  9) 
is  made.   However,  in  the  sensible-heat  model  visualized  in 
this  thesis,  a  linear  decrease  in  H„  with  pressure  was  as- 
sumed above  level  'x',  Fig.  10b.   Note  that  level  'x1  is  as- 
sumed to  occur  so  close  to  k  =  10  that  one  may  solve  (5-1) 
for 

Hr(k  =  9)  =  j  Hr  (Bulk).  (5-4) 

Solution  of  (5-4)  by  use  of  (5-1),  (5-2)  gives  the  result 

2  K  e 

CnV  T   +  — , -  2  K  Y 

T  =  -2-LJL Az Ic_  <  (5_5) 

X  2  K 

Az 

Equation  5-5  has  the  same  format  as  (5-3).   As  was  mentioned, 
Kaitala  [1974]  has  suggested  the  use  K*  =  2  x  105cm2sec 
in  (5-3),  and  this  would  lead  to  identical  values  as  the  use 
of  K*  =  105cm2sec_1  in  (5-5).   The  latter  K*  value  was  tested 
by  Warner  [1974]  who  felt,  with  Kaitala,  that  either  K*  value 
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was  only  an  approximate  choice  of  the  eddy  coefficient  re- 
quiring at  least  some  tuning,  with  a  theoretical  dependence 
upon  the  larger-scale  synoptic  parameters  over  the  ocean. 

The  physical  implications  of  condition  (5-4)  along  with 
the  linear  decrease  of  Hr(P)  with  pressure  to  zero  at  level 
k  =  8  is  depicted  in  Fig.  10b.   This  condition  is  compatible 
with  uniform  convergence  of  the  sensible-heat  flux  FL  of  the 
surface  layer  into  the  overlying  layer  (8,  10).   Having  uni- 
form convergence  from  level  x  (near  surface)  to  level  8  is 
equivalent  to  maintaining  a  "well  mixed"  boundary  layer  of 
thickness  from  the  surface  to  8. 

With  either  case,  Equation  5-3  or  5-5,  the  sensible  heat 
convergences  of  Fig.  10-a  and  10-b  were  considered  applicable 
at  level  k  =  9  and  were  considered  to  be  identical.   In  this 
thesis  the  value  of  2  K  in  (5-5)  was  based  upon  an  arbitrary 
choice  K*  =  105cm2sec    in  (5-3b),  rather  than  2  x  105cm2sec 

In  order  to  evaluate  H„ ,  it  was  necessary  to  derive  the 
surface  winds  V  by  use  of  (2-21)  at  all  gridpoints  of  the 
four  oceanic  meridians  under  study.   This  was  done  using 
Equations  2-16  and  2-20  with  contour  gradients  extracted  from 
the  1000  mb  field  except  in  the  latitude  band  ±5°.   Here,  the 
January  climatological  values  of  gradient  winds  were  taken 
from  the  Tropical  Wind  Atlas  of  Atkinson  [1970]. 

Besides  having  computed  the  applicable  surface  winds  V 
for  use  in  (5-2),  a  modification  to  C~  =  1.4  x  10  (after 
Weiler  and  Burling,  1967)  was  introduced,  as  compared  with 
CD  =  2 .  x  10    as  suggested  by  Kaitala  [1974]. 
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After  determination  of  T   in  Equation  5-5,  Hr  was  then 
computed  at  each  gridpoint  using  Equation  5-2.   The  values 
of  Hr  are  used  later  in  calculations  of  heat  budgets  in 
Sections  VI  and  VIII. 

C.   EVAPORATIVE  HEAT  TRANSPORT 

As  in  the  case  of  sensible  heat  convergence,  the  latent 
heat  transport  by  turbulence  was  assumed  to  be  subject  to 
uniform  heat  flux  convergence  throughout  the  layer  800-1000  mb 
(Fig.  lib).   The  total  amount  of  latent  heat  removed  by  evapo- 
ration from  the  ocean  surface  was  similar  to  that  given  by 
Kaitala  [1974]  using  bulk  transfer  theory.   While  a  turbulent 
model  similar  to  Fig.  lib  can  be  adopted,  it  is  necessary  to 
recognize  that  the  realized  latent  heat  may  occur  at  arbitrary 
levels  (for  example,  higher  than  k  =  8)  and  therefore  is  not 
introduced  at  level  k  =  9.   The  basis  of  the  turbulent  latent- 
heat  flux  model  adopted  here  is  shown  in  Fig.  lib,  where 

E(900)  =  WBulk)  (5-6) 

and  where  the  surface-layer  bulk  transfer  of  latent  heat  is 
given  by 

E  E  E(Bulk)=Lp1()  CD  Vs(q10-qx).        (5-7) 

From  the  planetary  boundary-layer  transfer  theory  analog  to 
Hr(900),  E(900)  may  be  written 

E(900)  =  Lp10  K  z9  -  z1Q  '  (5'8) 


In  (5-7),  and  (5-8),  L  is  the  latent  heat  of  vaporization 

io  cal  gm_1'  T10 


L  =  596.73  -  .601Tin  cal  gm"1,   Tin  in  °C. 
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Also  in  (5-8)  qg  and  q   are  mixing  ratios  at  the  level  k  =  9 
and  at  x,  the  top  of  the  constant  flux  layer.   The  factor  h 
which  appears  in  the  right  side  of  (5-6)  gives  the  type  of 
latent  heat  convergence  as  shown  in  Fig.  lib,  as  against 
that  of  Kaitala  (Fig.  11a).   Combining  Equations  5-6,  5-7, 
and  5-8  leads  to  the  solution  for  q  ,  namely 

CD  Vs  q10  +  2  h   q9 
q   =  — — - — ^-^ -  (5-9) 

X      C   V   +  2  K 
UD   s    L    ST 

with  Az  given  by  (5-3(c)).   As  in  the  case  of  T   in  Equation 
(5-5),  the  factor  of  2  k  of  (5-9)  was  taken  equal  to  Kaitala's 
K  of  Equation  5-3(b).   This  was  done  in  this  study  by  using 
a  value  of  K*  =  1.  x  105cm2sec    as  against  2  x  105cm2sec 
suggested  by  Kaitala.   Thus  amounts  computed  for  E  of  Equa- 
tion (5-7)  were  equivalent  to  amounts  computed  using  Kaitala's 
scheme  (Equation  5-8)  only  the  turbulent  convergence  of  E 
was  distributed  differently  in  the  planetary  boundary  layer. 
Thus  after  computing  q   (Equation  5-9),  E  was  calculated  at 

.X. 

each  gridpoint  except  those  over  ice  by  using  Equation  (5-7). 

D.   TURBULENT  HEAT  TRANSPORTS  OVER  AN  ICE -COVERED  OCEAN 

Any  gridpoint  whose  surface  temperature  was  less  than  -5°C 
was  considered  to  be  over  an  ice-covered  ocean  for  the  purposes 
of  calculating  HL  and  E.   In  this  study  there  were  four  such 
gridpoints.   In  the  handling  of  H„,  a  formula  adopted  from 
Kaitala  [1974]  was  first  utilized  to  describe  heat  conduction 
into  the  ice  (HICE) ,  from  the  warmer  water  beneath  the  ice 
surface: 

HICE  =  BBB*(T   -  T  )  (5-10) 

v  o    g 
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where  BBB  is  an  ice  conduction  coefficient  (697.83  ergs 

-  2      -  1  -  1 

cm   sec   (deg.  K)   ),  T   is  the  freezing  point  of  sea  water 
(271.3  deg.  K) ,  and  T   is  the  surface  temperature.   HICE  is 
then  utilized  in  the  following  surface  heat-balance  equation 
to  determine  the  summation  of  the  sensible  heat  flux  Hp  plus 
the  latent  heat  flux: 

Hp  +  E  =  HICE  +  QABG  -  F1Q*.  (5-11) 

A  Bowen  ratio  r  as  described  by  Kaitala  [1974]  was  then  cal- 
culated using  Equation  5-12,  where  the  right  side  is  an  em- 
pirical least-squares  best-fit  applicable  over  continental 
and  ice-covered  ocean  areas. 

H 
r(<{>)    e   -i-  =    9.6|Sin   <J>  |  2    -    7.93|Sin   c|>  |    +    2.0.     (5-12) 

Here  4>  is  the  latitude  of  the  gridpoint,  and  r  was  never  al- 
lowed to  exceed  1.33.   This  restraint  on  r  meant  that  E  was 
always  bounded  by 

E  <  .75  Hr  (5-13) 

(cf.  Kaitala,  1974).  From  the  definition  of  the  Bowen  ratio 
(r)  as  a  function  of  Sin  4> ,  and  using  the  following  equation 
for  the  combined  interface  fluxes  H   +  E  as 

Hr  +  E  =  Hr  (1  +  I  )  (5-14) 

one  obtains  simply 

Hr  =  ^  (Hr  +  E).  (5-15) 

Since  both  the  Bowen  ratio  (r)  and  Hp  +  E  are  known  at  each 
gridpoint,  Hp  is  easily  determined  from  (5-14).   Then  a  simple 
rearrangement  of  terms  in  the  Bowen  ratio  definition  permitted 
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the  computation  of  E.   This  is  given  by 

H 

E  ■  ?T57  •  (5"16) 

Equations  5-10  and  5-11  always  gave  negative  or  downward 
values  of  sensible  and  evaporative  heat  fluxes  at  the  four 
gridpoints  involving  ice-covered  oceans  since  QABG  -  F-  * 
was  dominantly  negative.   Further,  evaporation  never  exceeded 
75%  of  Hr  in  magnitude.   These  negative  values  appear  to  agree 
generally  with  those  arrived  at  in  the  work  of  Vowinckel  and 
Taylor  [1964]  which  dealt  in  part  with  evaporative  and  sen- 
sible heat  fluxes  over  an  ice-covered  Polar  Ocean  areas  in 
winter. 

E.   LARGE-SCALE  TURBULENT  HEAT  FLUX  ACROSS  THE  AIR- SEA  INTERFACE 

At  each  gridpoint,  it  is  possible  through  the  adaptation 
of  the  turbulent  boundary  layer  model  to  include  heating 
rates  across  the  air-sea  interface.   For  use  in  Section  VI, 
we  compute  the  values  of  H„  as  the  convergence  of  sensible- 
heating  rate  for  the  layer  (8,  10)  through  Equations  5-2,..., 
5-5.   Also  for  Section  VI,  the  combination  E  +  Hp  considered 
primarily  as  a  heat  loss  at  the  ocean  surface  is  computed  at 
each  gridpoint  by  use  of  Equations  5-7,  5-9.   For  the  over- 
ice  cases  Hp  +  E  is  determined  by  (5-11),  while  Hp  results 
from  (5-12)  and  (5-14).   In  any  case,  the  model  governing 
these  heating  rates,  Hp  and  E  has  been  extracted  from  the  FNWC 
operational  turbulent-boundary  layer  model  and  appended  to  the 
radiational  model  (with  the  2/3-CL  case)  for  operational  com- 
putations at  each  gridpoint  of  the  four  meridians  of  interest. 
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The  values  of  H„  turned  out  to  be  in  general  small  and 

negative  in  low  latitudes  reflecting  computed  negative  values 

of  T   -  T  .   However,  in  the  data-rich  mid-latitude  cross 
g    x  ' 

sections  where  both  strong  surface  winds  and  relatively  warm 
ocean  temperatures  coexist,  H„  is  positive  and  larger  and 
dominates  the  net  atmospheric  heating  in  layer  (8,  10).   Re- 
gardless of  the  inconsistency  in  the  sign  of  H„ ,  the  values 
of  E  +  Hr  seems  to  be  reasonable  heat-loss  rates  at  the  ocean 
surface  at  all  latitudes  insofar  as  comparison  with  other 
sources  of  data  (Budyko,  1956)  are  relevant. 

Finally  in  the  large-scale  sense  Hp  +  E  is  included  as  a 
(positive)  contribution  to  the  heat  budget  of  the  tropospheric 
column  even  though  it  is  not  known  precisely  at  what  level  the 
latent  heat  flux  E  will  be  realized. 


r\) 
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VI.   MERIDIONAL  CROSS-SECTIONAL  DEPICTION  OF  THE  HEAT  BUDGET 
SYSTEM  COMPUTATIONS 

A.   GENERAL 

The  general  design  of  this  section  is  to  utilize  all  com- 
putational subroutines  explained  in  Sections  III,  IV  and  V 
in  computations  of  a  single  time-step  of  the  heating  model 
developed  for  use  in  the  FNWC  prediction  model.   It  was  de- 
cided after  the  testing  of  F?  and  of  ALB  with  the  correspond- 
ing values  of  Raschke  et   al .     [1973]  that  only  the  two-thirds 
cloud  model  would  be  displayed  in  the  meridional  cross- 
sections  depicting  the  single  time-step  computations  of  the 
heating  package  of  the  FNWC  prediction  system  (Figs.  12, 
13(a),  13(b),  14(a),  14(b),  15(a),  15(b),  16(a),  and  16(b)). 
The  appropriate  radiative  calculations  with  two-thirds  CL 
were  thus  performed  at  each  gridpoint  of  the  four  meridians 
used  for  presentation  purposes  in  this  section.   However,  the 
same  type  computations  were  made  with  full-cloud  cover  so 
that  the  mean  radiative  balances  could  be  compared  layer-by- 
layer  and  at  the  levels  k  =  2  and  k  =  10.   Both  sets  of  com- 
putations were  zonally-averaged  and  used  to  get  mean  radia- 
tion meridional  cross -sections  displayed  in  Section  VII 
(Figs.  17,  18(a),  18(b),  19(a),  and  19(b)). 

In  this  study,  the  FNWC  gridpoint  processed  analyses  for 
0000GMT,  16  January  1974  were  used  at  the  three  Pacific  cross 
sections,  while  that  for  1200GMT,  16  January  1974,  was  used 
for  the  single  Atlantic  meridian.   In  any  case,  the  set  of 
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meridians  was  considered  as  representative  of  16  January 
data  when  converted  to  heating-rate  calculations. 

Figure  12  depicts  in  symbolic  language  the  key  to  the 
computational  entries  in  Figs.  13,  14,  15,  and  16.   In 
actuality,  this  key  list  represents  computations  made  at 
each  (I,  J)  gridpoint  sounding  from  data  in  the  form  of 
Table  1(b).   It  covers  those  computations  proceeding  from 
the  top  of  the  troposphere  to  the  ocean  surface,  including 
sensible  and  latent  heat  transfers  (Equations  5-2,  5-7). 
For  simplicity  and  ease  of  climatological  data  comparison, 
Figs.  13,  14,  15,  and  16  were  developed  by  interpolating 
gridpoint  results  to  whole  5°  latitude  increments.   The  in- 
terpolation routine  utilized  was  the  Lagrangian  cubic  inter- 
polation scheme.   It  has  the  form 

nm  -  0   (1-2)  (1-3)  (1-4)       (I-l)(I-3)(I-4)    f   n 
^UJ    ^1  (I-2)(I-3)  (1-4)    ^2  (2-1)  (2-3)  (2-4)    l°  ij 

+  0   (I-l)(I-2)(I-4)       (I-l)(I-2)(I-3) 
^3  (3-1) (3-2)  (3-4)    ^4  (4- 1) (4- 2) (4- 3) 

where  I  is  determined  from  an  equation  similar  to  (4-3(c)), 
Q(I)  is  the  desired  interpolated  value  of  the  heating  parameter 
at  a  non-integral  grid-value  I  that  corresponds  to  a  whole 
latitude  of  5°.   In  (6-1),  the  indices  1,  2,  3,  and  4  are 
for  any  four  increasing  successive  integral  values  of  I  along 
gridpoint  lines.   The  symbol  I  in  (6-1)  is  usually  applied 
between  points  2  and  3,  but  this  interpolation  scheme  may 
also  be  applied  with  good  accuracy  to  an  I-value  which  is 
non-centered  in  the  I-array  (1,  2,  3,  4).   In  order  to  com- 
plete the  interpolation  to  latitude,  the  array  is  applied 
sequentially  to  successive  sets  of  integral  I-values  along  each 
meridian.  q  y 


B.   EXPLANATION  OF  TERMS  IN  FIGURE  12 
1.   Cross-Section  at  Level  k  =  2 

It  should  be  explained  here  that  all  insolational 
parameters  enumerated  upon  in  Section  IV  dealt  with  the  speci- 
fic time  of  day  that  corresponded  to  the  hour  angle  h  at  the 
instantaneous  times  under  consideration.   The  resultant  in- 
cident solar  insolation  dealt  with  is  then  by  Equation  4-1 

F(2)  =  S(^-)~2  Cos  z. 
m 

In  order  to  avoid  reference  to  specific  map  times,  the  applic- 
able solar  hour-angles  were:   h  =  35°,  10°,  55°,  and  35°, 
respectively,  for  cross -sections  1,  2,  3,  4,  depicted  in 
Figs.  13,  14,  15,  and  16. 

Regarding  the  first  symbol  at  level  k  =  2  in  Fig.  12, 
note  QAVE  and  not  F(2)  is  displayed.   This  is  because  averages 
over  the  24-hour  day  (16  January  1974)  are  considered  more 
meaningful  for  the  data  day  under  consideration,  and  hence 
are  presented  in  this  form  in  Figs.  13,  14,  15,  and  16.   QAVE 
is  the  24  hour  average  of  F(2)  found  by  using  a  similar  form 


of  Equation  4-1  but  replacing  Cos  z  by  Cos  z  as  in  (6-2) 


QAVE  =  F(2)(  §§§— |-  )  (6-2) 


where 


Cos    z    =    [H   Sin    <j>    Sin    6    +    Cos    <J>   Cos    6    Sin  H]     (6-3) 


IT 

H(hour  angle  at  sunset)  =  arc  Cos[-  Tan  $   Tan  6] 

(f)  =  latitude 
6  =  solar  declination  for  16  January  1974 
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Cos  z   is  equal  to  the  24-hour  average  cosine  of  the  zenith 
angle.   The  24-hour  time-averaging  period  is  considered  con- 
sistent with  the  longwave  flux  divergences,  which  change 
slowly  with  the  time  of  day.   This  conversion  to  expected 
daily  averaged  heat- transfer  quantities  is  compatible  with 
the  determination  of  a  heat  budget  for  the  given  date  (e.g., 
16  January  1974) . 

In  keeping  with  the  concept  of  using  24-hour  averages, 
other  parameters  for  level  k  =  2  are  defined  as  follows: 


QREF  =  REF(t)(  — -|  )  (6-4) 

where 

REF(t)  =  F(2)  -  A26  -  A610  -  (IA10  +  IS10) .    (6-5) 

Here  REF(t)  is  the  instantaneous  solar-reflected  insolation 
at  the  gridpoint  as  given  in  4-30  and  QREF  is  its  24-hour 
average,  assuming  that  planetary  albedo  remains  constant  for 
the  24-hour  period,  even  though  it  was  calculated  for  the  in- 
dicated solar  time  at  each  meridian. 

The  same  principle  v/ill  be  used  with  regard  to  all 
other  solar  disposition  parameters  in  converting  from  time- 
dependent  values  at  solar  time  t  to  24-hour  averaged  values. 
Thus  in  the  following  Equations  6-7,  6-9  below  Q26(t)  be- 


comes Q26,  and  Q610(t)  becomes  Q610.   Superior  bars  (  )  will 
not  be  explicitly  shown  in  the  values  presented  in  the  cross- 
sections  key,  Fig.  12,  but  will  be  understood  in  what  follows 
Finally,  the  24-hour  average  "system"  balance  BALT  is  com- 
puted from 

BALT  =  QAVE  -  (QREF  +  F2*)  (6-6) 
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for  the  tropopause  level  k  =  2  at  the  indicated  latitude.   In 
computing  this  balance  the  value  of  F?*  was  introduced  from 
Equation  3-11.   Net  terrestrial  fluxes  MF*"  such  as  F~*  that 
appear  in  Fig.  12  were  considered  to  be  constant  throughout 
the  24-hour  period  which  would  be  approximately  true  if  the 
cloud-covers  were  representative  of  the  period. 
2 .   Cross -Section  in  Layer  (2,  6) 

Referring  to  Fig.  12  the  following  definitions  apply. 
In  this  layer   all  the  heat- transfers  shown  are  assumed  to 
be  radiative  only  in  our  version  of  the  FNWC  heating  model. 
The  averaged  radiative  cooling  (heating)  rate  is  given  by 


BAL26  =  Q26  -  F26  (6-7) 


where 


Q26  -  daily  average  absorption  in  layer  (2,6)  and  is 
defined  relative  to  A26(t)  by  a  transformation 

similar  to  Equation  6-2 

F26  =  terrestrial  cooling  rate  given  by  Equation 
3-12  (a). 

3 .   Cross-Section  in  Layer  (6,  10) 

When  considering  radiative  transfers  only  in  the  layer 

(6,  10)  the  24-hour  average  radiative  cooling  is  given  by 

BAL610  =  Q610  -  F610  (6-8) 

where  Q610  is  the  sum  of  Q68  and  Q810  in  Fig.  12.   BAL68  has 
been  taken  as  one-half  of  BAL610  .   However 

BAL810  =  Q810  -  F810  +  Hp  =  BAL68  +  H  .       (6-9) 

Here  Hf  is  the  sensible  heat  flux  transported  to  (8,  10)  by 
Equation  5-2.   Likewise,  F610  is  the  sum  of  F68  and  F810, 
terrestrial  flux  loss  values  for  the  respective  layers.   Note 
here  that  the  parameter  BAL610  remained  negative  at  all  gridpoints 
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of  Figures  13  through  16  when  considering  radiative  transfers 
only. 

4.   Cross-Section  at  Air-Sea  Interface  (k  =  10) 

The  heat  balance  at  the  earth's  surface  (BALB)  con- 
sisted of  the  following  equation 

BALB  =  QABG  -  F1Q*  -  (Hp  +  E) .  (6-10) 

Here  QABG  is  the  24-hour  average  insolation  absorbed  by  the 
surface  as 


QABG  =  QABG(t)  [Cos  z  /  Cos  z]  (6-11) 

(E  +  Hr)  is  the  sum  of  the  turbulent  heat  transfers 
which  were  computed  from  the  bulk  formulas  as  explained  in 
Section  V.B.  and  V.C. 

C.   MERIDIONAL  CROSS -SECTIONS  OF  THE  VERTICAL  HEAT  BUDGET 
FOR  16  JANUARY  19  74 


Figures  13,  14,  15,  and  16  as  explained  represent  the 
single-time  step  of  heating  computations  for  each  of  the  four 
meridians  used  in  this  study.   Note  that  the  four  Figures 
(13-16)  are  divided  into  a  and  b  parts  with  section  (a) 
showing  "tropical"  results  and  (b)  mid-to-high- latitude  re- 
sults.  While  the  results  depicted  in  these  cross -sections 
are  exhibited  as  representing  daily-averaged  values,  they  are 
actually  based  upon  heat-budget  parameters  at  the  specific 
map  time  0000GMT  and/or  1200GMT  on  16  January  1974.   Hence, 
if  the  results  are  to  be  used  as  part  of  the  "heat-package" 
subroutine  of  FNWC ,  all  solar-radiative  absorption  and-reflec 
tance  terms  must  be  recovered  as  functions  of  GMT  e.g., 


F(2,  t)  =  QAVE  *  [Cos  z  /  (Cos  z)] 
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Figure  14(a).   170°W  Longitudinal  cros 
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and 

QREF(t)  =  QREF  *  [Cos  z  /(Cos  z)  ] 

etc.   These  solar  disposition  terms  may  then  be  utilized  in 
connection  with  a  one-hour  time-step  application  of  the 
thermodynamic  equation  of  the  set  of  primitive  equations  used 
in  the  FNWC  prediction-process.   The  purpose  for  leaving  all 
results  of  Figures  13,...,  16  in  the  form  of  daily-averages 
was  primarily  so  that  these  figures  could  be  compared  with 
climatology  where  available  of  the  zonally-averaged  heat 
budget  of  the  ocean-atmosphere  system  for  January. 

Finally  while  only  the  (2/3) -CL  cloud  model  results  have 
been  depicted  in  Figs.  13,  14,  15,16.   An  analogous  set  of 
heat-budget  cross-sections  were  developed  for  the  full-CL 
model.   These  are  not  presented  here  because  of  the  sheer 
volume  of  the  results.   However,  certain  radiative  terms  in 
both  the  (2/3) -CL  and  full-CL  cross -section  sets  have  been 
extracted  from  both  sets  for  purposes  of  presenting  zonal 
comparisons  of  the  purely  radiative-heat-budget  of  the  FNWC 
data  as  it  was  applicable  to  16  January  1974.   The  compara- 
tive zonally-averaged  radiative  cross-sections  are  presented 
in  Figs.  18  and  19  of  Section  VII. 
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VII.   ZONALLY-AVERAGED  RADIATIQNAL  BALANCE  OF  THE  OCEAN- 
ATMOSPHERE  SYSTEM 

A.   GENERAL 

In  order  to  depict  the  results  of  the  purely  radiative 
contributions,  zonally-averaged  over  the  ocean-atmosphere 
system  for  this  study,  Figs.  18(a,  b)  and  19(a,  b)  are  pre- 
sented.  These  figures  present  respectively  the  zonally- 
averaged  radiative  cross -sect ions  for  both  the  2/3-CL  and 

s 
the  full-CL  cloud  amounts.   These  cross -sections  show  the 

results  after  zonal -averaging  over  the  four  meridians  con- 
sidered in  this  study.   The  results  are  then  displayed  in 
the  format  of  Fig.  17,  which  gives  the  key  to  the  line  items 
of  the  cross -sections .   Note  that  Fig.  17  simply  combines 
the  radiative-transfers  Q68,  Q810,  and  also  F68  with  F810 
of  Fig.  12  to  arrive  at  layer-average  radiational  warming 
rates  in  the  combined  layer  (6,  10).   Also  in  the  results 
of  this  section,  turbulent- transfer  heating  rates  into  the 
atmosphere  (and  corresponding  cooling  in  the  underlying  ocean) 
have  merely  been  omitted  from  the  results  of  Figs.  13,..., 16, 
in  passing  to  the  radiative  composites  presented  here. 

In  obtaining  a  simple  zonally-averaged  set  of  radiative 
heating-rates  at  each  latitude  cf>,  all  values  of  each  parameter 
listed  in  Fig.  17  at  latitude  $  in  the  range  20S,...,65N  (by 
5°  steps)  were  simply  averaged  over  the  meridians  providing 
the  averaged  value  at  <J> .   It  should  be  noted  that  at  <f>  =  65N, 
there  was  only  one  contribution  to  the  average.   In  the  Southern 
Hemisphere  region  20S,...,10S,  only  two  values  (from  125W  and 
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and  35W)  of  each  parameter  contributed  to  the  means  listed 
here  in  Figs.  18  and  19.   At  all  other  latitudes,  there  were 
in  general  four  values  of  each  radiative  parameter  entering 
into  the  computed  radiative  mean  cross -section.   As  a  re- 
sult, one  may  conclude  that  near  the  northern  and  southern 
boundaries  of  Figs.  18  and  19,  the  listed  values  of  the 
zonally-averaged  quantities  may  not  be  exact,  but  the  general 
zonal  trend  is  probably  reliable.   (Fig.  20  shows  this  to  be 
a  valid  result.) 

Figures  18,  19,  and  20  are  useful  for  the  discussion  of 
the  earth-atmosphere  system  radiative  balance,  and  in  the 
discussion  of  the  relative  merits  of  the  two-thirds  cloud- 
model  as  compared  with  the  full-CL  model.   The  full-cloud 
model  used  by  Warner  [1974],  is  also  used  in  comparative 
discussion  of  the  heat  budget  of  the  earth-atmosphere  system 
of  Section  VIII. 

B.   EARTH -ATMOSPHERE  SYSTEM  RADIATIONAL  BALANCE  RESULTS 

Figure  20  depicts  the  zonally-averaged  distribution  of 

three  parameters,  R  ,  R  ,  and  R  for  full  cloud  amounts  and 
r         '   s '   a' 

two-thirds  cloud  amounts,  respectively.  These  parameters 
were  defined  (after  Malkus ,  1962)  as  (1)  R  ,  the  mean  radia- 
tive energy  transfer  rate  across  the  top  of  the  ocean-tropo- 
spheric  system  (referred  to  as  BALT  at  k  =  2  in  Fig.  17;  (2) 

R  ,  the  mean  radiative  cooling  rate  in  the  troposphere,  i.e., 
a 


the  sum  of  BAL26  and  BAL610  in  Fig.  17;  and  (3)  R,  the  mean 
radiative  warming  (cooling)  rate  at  the  earth's  surface, 


referred  to  as  BAL10  in  Fig.  17.   The  relationship  between 
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10  20  ,  30 

Latitude   (°N 


Figure  20 


Radiational  balance  at  the  tropopause  (R  ),  at 


the  ocean  surface  (R),  and  in  the  tropospheric  column  (R  ). 

a 

Solid  lines  denote  computations  made  with  the  2/3-CL  (sub- 
scripted "1")  cloud  model  and  dashed  lines  correspond  to  the 
full-CL  model,  (subscripted  "2").   All  computations  are  for 
16  January  1974  on  ly(min)  x  102 . 
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these  three  radiative  parameters  is  as  follows: 

Rs  =  R  +  Ra.  (7-1) 

Table  V  lists  zonally-averaged  comparisons  of  R  ,  R  ,  and  R 

s   a 

as  a  function  of  latitude  <J>  for  this  study  for  16  January 
1974,  both  for  the  2/3-CL  and  the  full-CL  cases.   Also  in- 
cluded is  the  R  .    values  (signifying  satellite  climatology) 

after  Raschke  et   at.     [1973],  and  R7  and  R  -  values.   The 

i>  a  j 

latter  two  symbols  denote  December  climatological  values 
after  Berliand  [1956]  corresponding  to  R  and  R  respectively. 
"Hemispheric"  means  of  any  Q-value  are  then  by  the  following 
cosine-weighting  scheme: 

65 

•  E„  Q-  Cos  (J> . 

0        =   1  =  0   x 11'  f7-21 

^Wt.Avg.      65  *  l/  L) 

-En  Cos  d)  - 
1  =  0      Yi 

The  tabular  values  of  R  ,  R,  and  R   as  a  function  of 

s '   '      a 

latitude  and  for  the  two  cases,  2/3-CL,  and  full-CL,  have 

been  plotted  in  Fig.  20,  where  some  interesting  results  of 

the  radiational  model  become  evident. 

It  is  evident  from  Fig.  20  that  the  R  values  obtained 

to  s 

with  the  2/3-CL  cloud  model  exceed  those  based  upon  the  full- 
CL  model  up  to  <j>  =  40N.   Poleward  of  this  latitude  the  two 
R  -curve  cross  with  R  (CL)  becoming  slightly  larger  (although 
negative)  than  R  (2/3-CL) .   The  same  general  results  hold  for 

both  R.,  and  R~  and  for  R  -.    relative  to  R  ot    up  to  latitudes 

J.       z  a  1.  <-i  l 

40-45N,  where  also  cross-overs  in  the  graphical  depictions 
occur  in  both  cases.   These  results  are  related  to  the 
zonally-averaged  radiative  budgetary  cross -sect ions  of  Figs. 
18  and  19  as  is  explained  in  Section  VII  (C) . 
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Earth-Troposph 
Net  Radiation, 
(k=2) 

ere 

R 
s 

Surface  Net 

Net  Radiation 

(k=10) 

,  R 

Atmospheric-Net 

Radiation  Gain 

R 
a 

Lat 

R 

R 
si 

R 

S2 

Ri 

R2 

R3 

R 
ai 

R 

a2 

R 
a3 

20S 

.1800 

.0681 

.0216 

.1622 

.1123 

-.0942 

-.0913 

15 

.1700 

.0552 

.0052 

.1624 

.1125 

-.1073 

-.1075 

10 

.1600 

.0488  - 

.0039 

.1780 

.1310 

-.1292 

-.1347 

5 

.1430 

.0475 

.0046 

.1902 

.1537 

-.1426 

-.1492 

0 

.1150 

.0182  - 

.0270 

.1802 

.1443 

.1713 

-.1621 

-.1720 

-.0833 

5 

.0825 

-.0050  - 

.0478 

.1697 

.1387 

-.1746 

-.1864 

10 

.0550 

-.0140  - 

.0428 

.1645 

.1454 

.1481 

-.1787 

-.1882 

-.1157 

15 

.0150 

-.0489  - 

.0736 

.1374 

.1226 

-.1863 

-.1962 

20 

-.0125 

-.0813  - 

.1010 

.1123 

.1022 

.1088 

-.1937 

-.2031 

-.1551 

25 

-.0256 

-.0938  - 

.1058 

.0807 

.0738 

-.1745 

-.1796 

30 

-.0400 

-.0998  - 

.1041 

.0475 

.0429 

.0463 

-.1472 

-.1471 

-.1551 

35 

-.1050 

-.1497  - 

.1551 

.0152 

.0145 

-.1648 

-.1696 

40 

-.1375 

-.1689  - 

.1706 

-.0092 

-.0072  - 

.0116 

-.1597 

-.1634 

-.1597 

45 

-.1625 

-.1884  - 

.1835 

-.0094 

-.0068 

-.1791 

-.1768 

50 

-.1750 

-.2062  - 

.1969 

-.0275 

-.0200  - 

.0417 

-.1788 

-.1769 

-.1828 

55 

-.2025 

-.1798  - 

.1671 

-.0409 

-.-0300 

-.1388 

-.1371 

60 

-.2300 

-.1932  - 

.1759 

-.0868 

-.0683  - 

.0509 

-.1065 

-.1077 

-.2037 

65N 

-.2400 

-.2265  - 

.2265 

-.0868 

-.0868 

-.1397 

-.1397 

Wt. 
Avg 
(N. 
Hem) 

-.0500 

-.1000  - 

.1139 

.0668 

.0581 

.0699 

-.1668 

-.1720 

-.1437 

Table  V.   Zonally-averaged  components  R,  R  ,  R  of  the  earth-atmosphere 

S     3. 

system  as  a  function  of  latitude.   The  averaged  values  as  found  in  this 
study  using  subscript  "1"  for  the  2/3-CL  case  and  the  subscript  "2"  for 
the  full-CL  case.   Corresponding  values  of  R  reported  by  Raschke  et  al. 
(1973)  are  denoted  by  subscript  "4"  and  those  of  R  and  R  reported  after 

3 

Berliand  (1956)  for  December  climatology  denoted  by  subscript  "3". 
Weighted  averages  over  the  latitude  range  0-65N  are  shown  in  the  last 
line  of  the  table.   (All  values  in  ly  min-1). 
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C.   EXPLANATION  OF  THE  MEAN  ZONAL  RADIATIVE-BUDGET  DISTRIBUTIONS 

It  should  be  observed  that  in  Table  V,  the  R  values  ob- 

'      s 

tained  with  the  2/3-CL  model  are  greater  than  those  resulting 
from  the  full-CL  model  throughout  the  tropics  and  up  to  lati- 
tude 40N  (approximately).   Poleward  of  latitude  40N,  the 
values  of  R  (2/3-CL)  become  smaller  than  those  for  full-CL. 
These  comparisons  are  illustrated  clearly  in  Fig.  20  by  the 
zonal  depictions  of  the  R  values  calculated  by  the  two 
cloud-model  cases.   Examination  of  Table  V  for  comparative 
values  of  the  radiative  balances  at  the  surface  shows  the 
same  general  effect  at  the  ocean  surface,  with  the  cross-over 
between  2/3-CL  and  full-CL  cases  (R,  and  R2  curves  in  Fig.  20) 
occurring  near  latitude  40N.   Finally,  the  same  general  effect 
is  evident  in  the  comparative  values  of  R  (2/3-CL)  compared 

to  R  (CL) .   That  these  three  comparisons  should  follow  each 
a 

other  systematically  is  perhaps  to  be  expected  in  view  of 
the  connecting  relationship 

R   =  R  +  R  . 
s        a 

Of  the  changes  resulting  from  substituting  (2/3-CL)  for  the 
full-CL  cases,  the  largest  gain  in  low  and  mid- latitudes  oc- 
curs in  the  comparative  R  -curve  for  the  (2/3-CL)  case.   The 

smallest  gain  occurs  for  the  computed  R  -changes,  and  an  in- 

a 

termediate  gain  is  registered  for  the  R(2/3-CL) ,  the  net 
radiative  flux  at  the  ocean  surface. 

The  same  results  for  R  -changes  may  be  seen  by  comparing 
Fig.  18  for  the  2/3-CL  case  with  Fig.  19  for  the  full-CL  case. 
In  such  a  comparison  it  is  seen  that  the  systematic  reduction 
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of  CL  to  2/3-CL  increases  the  net  useful  solar  insolation 
at  k  =  2,  (Figs.  17,  18,  and  19)  by  reducing  QREF.   However, 
the  reduction  achieved  in  QREF(2/3-CL)  is  quite  small  at 
latitudes  <$>    >  45N.   On  the  other  hand  with  the  2/3-CL  case, 
the  F-*  values  (net  terrestrial  loss  to  space  at  k  =  2)  is 
increased  by  a  small  but  consistent  amount  which  at  higher 
latitudes  <$>    >    45N,  exceeds  the  very  small  net  solar  insola- 
tion gain  (k  =  2) . 

It  therefore  turns  out  that  in  the  consideration  of 
Figs.  18  and  19  of  the  atmospheric-layer  radiative  transac- 
tions (both  solar  and  terrestrial)  that  the  net  radiation 
R  penetrating  to  the  surface  is  influenced  primarily  by  QREF, 
at  best  in  lower  latitudes  and  up  to  cf>  =    40N.   Here  -A (QREF) 
exceeds  A  F  *  as  a  result  of  the  change  in  cloud  cover  and 

thus  accounts  for  greater  R  -values.   The  same  effects  ac- 

°        s 

count  for  R-changes  at  the  surface,  while  at  <f>  >  45N  the 

cross-over  effect  holds  both  at  k  =  2,  and  also  at  k  =  10. 

The  curves  of  R  versus  cf>  in  Fig.  20,  in  effect  measure 

a 

the  atmospheric  participation  in  the  radiative  absorption  and 
emission  processes,  follow  the  lead  of  -A(QREF)  in  lower 
latitudes  and  of  F  *  in  higher  latitudes,  which  act  as  the 
forcing  functions  in  the  radiative  balance  of  the  system. 

Recent  satellite  results  especially  of  Von  der  Haar  and 
Hanson  [1969]  as  well  as  of  Raschke  et   at.     [1973] ,  also  draw 
attention  to  the  reduced  values  of  global  albedo,  compared 
to  previous  heat  budget  estimates  (e.g.,  London,  1957)  in 
tropical  and  subtropical  latitudes.   This  also  shows  up  in 
Table  V  in  the  comparative  values  of  Raschke' s  distribution 
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of  R  with  the  "best"  R   distribution  computed  here  for  the 

2/3-CL  model.   It  should  be  noted  from  Table  V  that  the  use 

of  the  2/3-CL  model  leads  to  R  -,  -values  which  are  in  closer 

si 

agreement  with  R  .    up  to  latitude  40N,  than  that  generated 

by  the  full-CL  model.   At  high  latitudes  in  Table  V  the 

comparisons  between  the  results  of  these  two  cloud  models 

are  inconclusive  in  their  relationship  to  Raschke's  results. 

However,  his  results  in  some  of  these  areas  are  over  ice, 

which  enhances  his  effective  QREF  as  compared  with  either 

model  used  here. 

Finally  note  that  comparative  values  after  Berliand 

[1956]  are  also  included  in  Table  V  for  R  and  R   at  10°  lati- 
L     J  a 

tude  intervals  based  on  December  global  climatology.   The 
most  noteworthy  comparisons  to  be  made  here  is  the  transi- 
tion from  positive  to  negative  surface  net  flux  R,  which  oc- 
curs between  35-40N,  in  reasonable  agreement  with  the  results 
interpolated  either  from  the  R. - ,  R~ -distribution  of  Table  V. 
In  general,  the  global  results  of  Berliand  reflect  a  greater 
continental  trend  with  latitude  than  is  found  in  the 
January  1974  study  over  the  oceans. 
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VIII.   ZONALLY-AVERAGED  TROPOSPHERIC  AND  OCEANIC  HEAT  BUDGETS 
FOR  16  JANUARY  19  74 

A.   THE  TROPOSPHERIC  HEAT  BUDGET  OVER  ICE -FREE  OCEANS 

By  averaging  over  the  tropospheric  column  for  the  four 

meridians  displayed  in  Figs.  13,..., 16,  latitude  by  latitude, 

the  zonally-averaged  tropospheric  balance  results  (Table  VI 

and  Fig.  21).   Values  presented  in  Table  VI  and  Fig.  21  are 

R  ,  the  tropospheric  radiative  net  cooling  defined  in  Section 
a 

VIII. B.,  and  E  +  Hr  as  formulated  in  Equations  5-2,  5-7.   The 
latter  parameter  represents  the  combined  turbulent  heat-trans- 
fer Hy    +  E  across  the  ocean-air  interface  (see  for  example, 
Figs.  13,..., 16).   The  heat-budget  equation  for  the  atmosphere 
column  at  latitude  <f>  then  becomes,  with  the  notation  of 
Malkus  [1962]: 

Qva  +  Sa  =  Ra  +  <E  +  Hr^ ■  (8-1) 

Here  S   is  the  in-situ   storage-heating  rate  of  the  column, 
a 

and  Q    is  the  required  flux-divergence  of  heat  necessary  to 

bring  about  the  "balance"  in  the  column,  for  the  observed 

values  of  R  E  +  H„  and  S  . 
a      r      a 

Table  VI  lists  the  zonally-averaged  values  of  R  ,  E  +  Hp 
and  the  resulting  values  of  Q   +  S  ,  which  are  computed  from 
(8-1).   As  before  the  subscript  notation  "1"  and  "2"  denote 
the  2/3-CL  and  full-CL  computations  respectively.   Figure 
21  depicts  each  function  R  ,  E  +  H„  and  Q   +  S  ,  computed 
by  the  two  cloud-models. 

As  noted  before  in  Section  VII. C,  R  ,    tends  to  exceed 

'   al 

R  -  over  most  of  the  tropical  Southern  Hemisphere  and  in  the 
a  z 
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r 

equirement 

s 
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a 
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Q  +S 
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R 
ai 

R 
a2 

R 

a3 

(E+Hr) 
*  l 

(E+Hr) 

1  2 

(E+Hr) 

1   3 

(Q  +s  ) 
va  a 

(Q  +s  ) 
i   va  a  2 

(Q  +S  ) 
va  a  3 

20S 

-.0942 

-.0913 

.1402 

.1402 

.0461 

.0490 

15 

-.1073 

-.1075 

.1718 

.1718 

.0644 

.0642 

10 

-.1292 

-.1347 

.1689 

.1689 

.0397 

.0342 

5 

-.1426 

0.1492 

.1197 

.1197 

-.0229 

-.0294 

0 

-.1621 

-.1720 

-.0833 

.0839 

.0839 

.0764 

-.0782 

-.0882 

-.0069 

5 

-.1746 

-.1864 

.0791 

.0791 

-.0956 

-.1073 

10 

-.1787 

-.1882 

-.1157 

.1374 

.1374 

.1018 

-.0413 

-.0508 

-.0139 

15 

-.1863 

-.1962 

.1611 

.1611 

-.0252 

-.0351 

20 

-.1937 

-.2031 

-.1551 

.1430 

.1430 

.1342 

-.0507 

-.0601 

-.0209 

25 

-.1745 

-.1796 

.1257 

.1257 

-.0488 

-.0538 

30 

-.1472 

-.1471 

-.1551 

.1998 

.1998 

.1250 

.0526 

.0528 

-.0301 

35 

-.1648 

-.1696 

.2341 

.2341 

.0692 

.0644 

40 

-.1597 

-.1634 

-.1597 

.2971 

.2971 

.1134 

.1374 

.1337 

-.0463 

45 

-.1791 

-.1768 

.2172 

.2172 

.0382 

.0404 

50 

-.1788 

-.1769 

-.1828 

.1901 

.1901 

.1458 

.0113 

.0132 

-.0370 

55 

-.1388 

-.1371 

.1301 

.1301 

-.0087 

-.0070 

60 

-.1065 

-.1077 

-.2037 

.1156 

.1253 

.1944 

.0091 

.0172 

-.0093 

65N 

-.1397 

-.1397 

-.0721 

-.0721 

-.2118 

-.2118 

Wt. 

Avg 
(N. 

-.1668 

-.1720 

-.1437 

.1508 

.1512 

.1208 

-.0160 

-.0208 

-.0229 

Hem) 

Table  VI.   Zonally-averaged  components  R  ,  (E+Hr)  and  (Q  +S  )  of  the 

a      I        va  a 

tropospheric  heat  budget.   The  averaged  values  as  found  in  this  study 
using  subscript  "1"  for  the  2/3-CL  case  and  the  subscript  "2"  for  the 
full-CL  case.   Corresponding  values  reported  after  Berliand  (1956)  for 
December  climatology  denoted  by  subscript  "3".   Weighted  averages  over 
the  latitude  range  0-65N  are  shown  in  the  last  line  of  the  table.   (All 
values  are  in  ly  min-1). 
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Figure  21.   Distribution  of  R   E+Hr  and  of  Q  ,  +So  as  a 

a      i  va   a 

function  of  latitude  for  the  two  cloud-model  cases. 

(solid  curves  —  2/3-CL  case;  dashed  curves  —  full-CL  case) 
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Northern  Hemisphere  to  close  to  40N,  although  these  differences 

are  quite  small.   The  distribution  of  E  +  L  as  a  function  of 

latitude  is  independent  of  the  cloud-model  except  at  high 

latitudes  over  a  frozen  sea-surface  which  occurs  near  <f>  >  60N. 

Figure  21  shows  the  identical  sets  of  E  +  H„  up  to  latitudes 

55N.   The  curves  of  Q   +  S   are  also  shown  in  Fig.  21,  as 

xva    a  6     * 

the  residual  of  the  right  sides  of  Equation  8-1  for  both  the 

2/3-CL  and  full-CL  cases,  respectively.   Note  that  (Q   +  S  ) - 

'     r       3  vxva    a.J  1 

slightly  exceeds  (Q,ro  +  S  )9  over  most  of  the  latitude  range 

va    a  z 

to  25N  as  a  result  of  the  corresponding  of  excess  R  ,  over 

R  9  for  the  same  latitude  range.   Note  also  that  both 

Q   +  S   distributions  reach  peak  values  at  <t>  =  40N,  and 

that  the  values  Q    +  S   are  generally  positive  over  30-50N. 

xva    a     &       7  r 

This  latter  feature  follows  from  (8-1),  since  Q    +  S 

*■    J  '       xva    a 

must  follow  the  trend  of  E  +  Hr  where  the  latitudinal  varia- 
tion of  the  latter  function  peaks  sharply,  relative  to  more 


or  less  uniform  behavior  of  R  (<j>) 

a 


The  hemispheric  mean  value  of  (Q    +  S  )  =  -.0160  Zy   min 
F  ^xva    a.J  J 


-  i 


per  cm   column.   In  this  study,  Q  .  has  not  been  computed,  but 

va 


if  the  (Q    +  S  )  is  attributed  to  mean  storage-cooling  S 
va    a  a 

of  the  troposphere  alone,  the  daily  storage  rate  corresponds 
to  a  temperature-change  rate  given  by 


w          (R   +  E  +  H  ) 
<  ll  ^a=  4'X     AP  mb ~   x  1440  min/da^    ^8-2^ 

with  AP  mb  =  800  mb  in  the  troposphere.   This  resultant  mean 
in-situ   hemisphere  cooling  rate  (8-2)  over  the  tropospheric 
depth  of  800  mb  is  then 
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fl)a  -  -.12»C(dayr' 

averaged  over  the  cm2  tropospheric  column.   This  seems  to  be 
a  reasonable  midwinter  global  value. 

The  variation  in  the  computation  of  E  +  H„  in  the  2/3-CL 
and  the  full-CL  cases  which  occurs  poleward  of  latitude  55N 
is  explained  in  Section  VIII. C. 

B.   THE  ZONALLY -AVERAGED  HEAT  BUDGET  OF  THE  ICE -FREE  OCEAN 

For  the  present,  neglect  the  minor  variations  attributed 
to  water-to-ice  conduction  rates  which  occur  in  latitudes 
poleward  of  55N  (see  for  example  H.    of  Equation  5-10). 

1  VV 

Neglecting  this  ice-conduction  term,  the  zonally-averaged 
heat-budget  of  the  ocean  water-mass  obeys  the  thermodynamic 
relationship 

Qvo  +  S  =  R  +  [-(E  +  Hr)]  (8-3) 

for  both  the  2/3-CL  and  full-CL  cases.   Note  that  Equation 
8-3  is  the  analog  of  (8-1)  for  the  atmosphere,  with  -  (E+H„) 
representing  generally  the  cooling  rate  for  the  water  mass. 

The  terms  of  the  right  side  of  (8-3)  have  been  computed 
for  each  meridian  and  averaged  by  latitude,  with  the  results 
shown  in  Table  VII.   Then  Q   +  S  is  computed  as  the  residual 
of  the  sum  of  the  two  functions  on  the  right  side  of  (8-3), 
and  shown  also  in  columns  1  and  2  of  Table  VII.   The  graphi- 
cal distribution  of  the  results  listed  in  this  table  are 
shown  in  Fig.  22  both  for  the  2/3-CL  (subscript  "1")  and 
the  full-CL  (subscript  "2")  cases. 

As  before  the  primary  ocean-mass  heating- function  is  R, 
the  radiative  heating  rate  at  the  ocean-air  interface.   For 
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Surface  Net 

Radiation,  R 

(k=10) 

Turbulent  transport 
heating  rate  loss 

Surface  Balance 
Requirements 

Lat 

Ri 

R2 

R3 

-(E+Hr)!-(E+Hr)2 

~(E+Hr)3 

((Vo+S) 

>%o+S) 

;(Qvo+S)3 

20S 

.1622 

.1123 

-.1402  -.1402 

.0220 

-.0280 

15 

.1624 

.1125 

-.1718  -.1718 

-.0094 

-.0592 

10 

.1780 

.1310 

-.1689  -.1689 

.0091 

-.0379 

5 

.1902 

.1537 

-.1197  -.1197 

.0705 

.0340 

0 

.1802 

.1443 

.1713 

-.0839  -.0839 

-.0764 

.0964 

.0604 

.0949 

5 

.1697 

.1387 

-.0791  -.0791 

.0906 

.0596 

10 

.1645 

.1454 

.1481 

-.1374  -.1374 

-.1018 

.0271 

.0080 

.0463 

15 

.1374 

.1226 

-.1611  -.1611 

-.0238 

-.0385 

20 

.1123 

.1022 

.1088 

-.1430  -.1430 

-.1342 

-.0306 

-.0408 

-.0254 

25 

.0807 

.0738 

-.1257  -.1257 

-.0450 

-.0520 

30 

.0475 

.0429 

.0463 

-.1998  -.1998 

-.1250 

-.1524 

-.1570 

-.0787 

35 

.0152 

.0145 

-.2341  -.2341 

-.2188 

-.2195 

40 

-.0092 

-.0072  - 

.0116 

-.2971  -.2971 

-.1134 

-.3063 

-.3043 

-.1250 

45 

-.0094 

-.0068 

-.2172  -.2172 

-.2266 

-.2240 

50 

-.0275 

-.0200  - 

.0417 

-.1901  -.1901 

-.1458 

-.2176 

-.2101 

-.1875 

55 

-.0409 

-.0300 

-.1301  -.1301 

-.1710 

-.1601 

60 

-.0868 

-.0683  - 

.0509 

-.1156  -.1253 

-.1944 

-.2024 

-.1933 

-.2453 

65N 

-.0868 

-.0868 

.0721   .0721 

-.0147 

-.0147 

Wt. 
Avg 
(N. 

Hem) 

.0668 

.0581 

.0699 

-.1508  -.1512 

-.1208 

-.0840 

-.0931 

-.0509 

Table  VII.   Zonally-averaged  components  R,  -(E+Hf)  and  (Q  +S)  of  the 
surface  heat  budget.   The  averaged  values  as  found  in  this  study  using 
subscript  "1"  for  the  2/3-CL  case  and  the  subscript  "2"  for  the  full-CL 
case.   Corresponding  climatology  values  for  the  month  of  December  re- 
ported after  Berliand  (1956)  are  denoted  by  subscript  "3".   Weighted 
averages  are  listed  in  the  final  line  for  the  latitude  range  available 
in  the  N.  Hemisphere.   (All  values  in  ly  min   .) 


129 


(Qvo+Sdl 


(Qvo+s, 


(E  +  Hr)i&2 


20  10 

°S 


10  20    .         30 

Latitude   (°N 


40 


50 


60 


Figure  22.   Distribution  of  R,  -  (E+H  )  and  of  (Qvo+SQ) 
as  a  function  of  latitude  for  the  two  cloud-model  cases, 
(solid  curves  —  2/3-CL  case;  dashed  curves  --  full-CL  case) 
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the  2/3-CL  case,  the  corresponding  R1  is  considerably  larger 
in  low-  and  mid- latitudes  than  R?  ,  for  the  full-CL  case. 
There  is  virtually  no  change  in  the  - (E  +  H  )  profiles  over 
latitude  (except  poleward  of  60N) .   This  means  that  the  in- 
ternal heating  rate  of  the  ocean-mass,  which  is  given  by  the 
sum  of  the  two  terms  on  the  right  side  of  (8-3),  is  generally 
greater  in  the  2/3-CL  case  than  in  the  full-CL  case.   Here 
Q    is  the  required  mean  oceanic  heat-flux  divergence  at 
latitude  cj>.   Figure  22  supports  the  conclusions  of  Von  der 
Haar  and  Oort  [1973]  which  requires  increased  lateral  trans- 
port-divergences Q    in  the  oceans,  based  upon  the  recently 
observed  satellite-values  of  the  global  albedo.   The  values 

ascribed  here  to  (Q    +  S  ) ..  seem  to  corroborate  the  satel- 

^xvo    oJ  1 

lite  findings. 


The  "hemispheric"  mean  value  of  TO +  S  )   computed 
r  ^xvo     oJ  i  e 

between  0  -  65N  is 


-  i 


(Q   +  S)-  =  -.0859  ly   min 
xvo     J 1  } 

which  corresponds  to  negligible  cooling  rate  in  the  water-mass 
column,  although  it  is  five  times  as  great  as  the  correspond- 
ing atmospheric  effect. 

C.   OCEAN-TROPOSPHERIC  HEAT-BUDGETS  FOR  ICE-COVERED  OCEANS 
1 .   Computation  of  (E  +  H?)  Over  Ice-Covered  Oceans 

The  determination  of  E  +  H„  over  ice-covered  ocean 
was  computed  (following  Equation  A-36  of  Kaitala,  1974)) 
with  the  following  interface  conditions  used  in  Section  V: 

E  +  HP  =  R  +  H.  (5-11) 

r        ice  v    ' 

H.    =  BBB(T  -T  ) .  (5-10) 

ice       *-  o   g' 
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In  latitudes  where  ice-covered  oceans  occur,  R  is  negative  in 

January,  and  H.    is  usuallv  positive  but  of  an  order  of 
J  ice  '  r 

magnitude  smaller  than  the  magnitude  of  R.   Hence  the  (E  + 

H„)  values  computed  at  both  <j>  =  60N  and  <j>  =  65N  from  the 

170W  meridian  are  both  negative  (here  the  Bering  Sea  is 

frozen  on  January  16,  1974). 

Solution  for  H.    and  Hp  +  E  at  d>  =  60  and  65  N  then  gave 
ice      r  b 

the  following  results 


TABLE  VIII.   Surface  heat- transports  over  ice  at  170W 
Lat  Hice  iy   min'1  E  +  Hp  £y  min"1 


2/3-CL   FULL-CL  2/3-CL   FULL-CL 

<J>  =  65N      .0147     .0147  -.0721   -.0721 

4  =  60N      .0041     .0070  -.0683   -.0513 

The  values  of  E  +  H„  of  the  foregoing  display  are  used  in  the 
computation  of  zonally-averaged  heat  balance,  both  in  the 
troposphere,  and  the  ocean  as  reproduced  in  Tables  VI  and 
VII,  respectively. 

The  reason  for  the  difference  in  computed  values  of  E  +  Hp 
associated  with  the  2/3-CL  and  full-CL  cases  at  latitudes  $ 
>    60N  is  associated  primarily  with  Equation  5-11  and  the 
presence  of  ice.   At  <J>  >  60N,  the  surface  radiative  quantities 
satisfy  the  inequality 

|Rj  <  |  R-  |  ,  when  CL  >  0 

although  both  are  negative.  Furthermore,  H.  in  either  case 
is  an  order  of  magnitude  smaller  than  R^  or  R, ,  respectively. 
Hence,  when  CL  >  0 
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|E  ♦  Hr|2  <  |E  ♦  Hplj 

and  there  must  be  a  smaller  turbulent  flux  from  air  to  ice 

in  the  full-CL  case.   The  results  of  Table  VIII  at  latitude 

60N  exemplify  this  feature. 

2 .   Budgetary  Equations  in  Both  Media  with  Ice-Covered 
Ocean 

Finally,  H.    represents  the  result  of  a  molecular 
' '   ice    r 

conduction  transfer  rate,  which  leaves  the  ocean  and  enters 
the  atmosphere  over  the  ice  by  a  process  not  included  in  the 
turbulent  exchange  rate  E  +  Hr .   Hence  at  ice-covered  grid- 
points  the  heat  budget  equations  of  the  troposphere  and 
water-body  become: 

Air-mass  over  ice: 

Q    +  S   =  R   +  (E  +  Hr)  +  H.  (8-4) 

xva    a    a    ^     r^     ice  v   J 

Water-body  under  ice: 

Q    +  S   =  R  -  (E  +  Hr)  -  H.   .  (8-5) 

xvo    o        y  Y J  ice  v    J 

The  additive  terms  on  the  right  side  of  8-4,  8-5  did  not 

modify  the  weighted-mean  heating  rates  of  the  air-column  or 

the  water  body  significantly  over  the  values  resulting  from 

Figs.  21  and  22. 
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IX.   CONCLUSIONS 

The  present  work  has  been  primarily  involved  with  the 
development  of  the  radiational  two-layer  cloud  model  adapted 
from  Warner's  earlier  thesis  [1974].   The  data  employed  was 
oceanic  data  from  the  FNWC  gridpoint  analyses  for  16  January 
1974.   Here,  it  was  found  that  the  net  radiative  flux  avail- 
able for  the  earth- troposphere  system  was  best  verified 
against  January  satellite  climatology  (Raschke  et    at.  ,  1973) 
when  the  Smagorinsky  cloud-amounts  were  arbitrarily  reduced 
by  a  factor  of  one-third.   It  is  recognized  that  there  were 
some  minor  differences  in  the  comparisons  which  could  be  at- 
tributed to  the  one  day  space-time  average  as  opposed  to  the 
average  of  several  weeks  for  satellite  data,  and  years  inso- 
far as  surface  climatological  data  was  concerned. 

Reduction  of  Smagorinsky  cloud  amounts  by  one  third  was 
found  to  provide  more  realistic  global  albedo  estimates  than 
resulted  using  unreduced  cloud-amounts.   In  the  comparison 
with  satellite  data,  the  radiative  model  with  a  reduced  CL 
gave  slightly  higher  albedo  estimates  in  the  tropics  and 
mid- latitudes  and  slightly  lower  albedo  in  the  high-latitudes. 
The  radiative  model  still  lacks  an  effective  surface  albedo 
model  over  ice. 

With  the  reduced  cloud  amounts  greater  surface  net  heating 
rates  in  the  low  to  mid  latitudes  were  computed  than  when 
using  the  original  Smagorinsky  cloud  amounts.   This  latter 
finding  is  in  accord  with  recent  heat  balance  studies  using 
satellite   data  which  have  found  an  increase  in  surface 
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heating  particularly  in  the  tropics.   It  is  felt  that  further 
investigation  of  the  tropical  radiative  balances  are  neces- 
sary.  Thought  should  also  be  directed  towards  the  geometry 
of  solar  radiation  reflection  problem  among  cumulus  subgrid 
cells  in  the  tropics  as  contrasted  with  the  larger-scale  two- 
cloud  layer  model  used  in  this  study,  which  proved  fairly 
reliable  in  middle  and  high  latitudes.   Moreover,  specifica- 
tion of  the  primarily  cumuliform  cloud  distributions  in  the 
tropics  is  evidently  closely  related  to  the  subgrid  scale 
convection  model  rather  than  to  the  large-scale  distribution 
of  relative  humidity  (of  Smagorinsky) .   It  is  recommended  that 
further  testing  of  the  radiative  model  for  the  different  sea- 
sons of  the  year  continue,  particularly  with  the  one-third 
or  similar  reduction  in  the  Smagorinsky  cloud-amounts. 

Further  differences  from  climatological  data  might  be 
attributed  to  the  treatment  of  cloud  reflectivity  and  absorp- 
tivity for  solar  wavelengths,  as  well  as  the  parameterization 
of  the  sea-air  exchanges  of  sensible  and  latent  heat.   In  the 
latter  case  a  new  interpretation  of  the  simple  bulk  formulas 
and  of  a  redefined  distribution  of  the  turbulent  sensible  and 
latent  fluxes  in  the  planetary  boundary  layer  parameterization 
model  seemed  to  give  a  new,  conceptually  satisfying  view  of 
the  planetary  boundary  layer  model.   Additional  testing  with 
the  redefined  planetary  boundary  layer  or  any  new  interpre- 
tation of  such,  might  also  prove  productive  in  providing  bet- 
ter comparisons  with  recent  investigations  such  as  the  GATE 
and  BOMEX  experiments  involving  surface  heat  budget  calculations 
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APPENDIX  A:   EMISSIVITY  e   FORMULAS  AS  A  FUNCTION  OF  WATER 
VAPOR  AND  Co!JCABSORBER  MASSES  ALONG  PATH  OF 
INTEGRATION. 


The  general  form  of  the  calculation  required  for  computing 
net  IR  flux  is: 

Bio 

F10*  "  B10  -  B£0    S,c(U>C>T)  dB  (A"1) 

where  B1Q  =  Stephen-Boltzman  blackbody  flux  at  T, Q  and  e 
is  the  emissivity  of  the  atmosphere  as  a  function  of  both 
water  vapor  and  CO^  absorber  masses  along  the  path  of 
integration.   The  hatched  areas  in  Section  III  figures  are 

represented  by    J         e   dB  which  is  the  key  integral  to  be 

B=0   wc 

computed.   In  order  to  perform  this  calculation ,  Sasamori 
[1968]  proposed  emissivity  formulas  of  the  following  form: 


e   (6,8)  =  (.240  log,n[U(6,8)]  +  .622) 
wc  10     ' 


+  .07262((1.  -  .62556[U(6,8)  +  .0286]'26) 


[loglQ  C(6,8)  +  1.064]}  (A-2) 


for  the  layer  6-8  relative  to  the  level  6  as  reference  level 
^or  temperatures  T  j>  210  deg  k  where  e   is  very  nearly 
temperature  independent.   Also  note  that  U  is  the  water 
vapor  absorber  mass  given  by  (2-7)  measured  from  the 
reference  level  10,  6,  or  2  and  C  is  the  correspondingly 
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carbon  dioxide  absorber  mass.   Water  vapor  emissivity  alone, 
contribution  is  represented  by  the  quantity  within  the  brace 
of  Eq .  A-2.   The  second  brace  indicates  the  additional  contri- 
bution made  by  the  part  of  the  COp  emissivity  transmitted 
through  the  water  vapor  "overlap"  near  15um,  i.e.,  the  part 
not  already  counted  in  this  spectral  region  by  the  water 
vapor  integration. 

Equation  A-2  was  used  to  compute  atmospheric  emissivities 

over  all  the  layers  (8,  10),  (6,  10),  (4,  10) (1,  10) 

for  which  the  running  mass  values  U(8,  10),  C(8,  10),  etc. 
were  ready  to  be  implied  from  the  radiative  soundings.   To 

complete  the  F   4-  integration  given  by  (A-l)  from  level  k  =  1 
wc 

to  the  origin  (B  =  0,  e    =  1.0),  on  the  Yamato  Radiation 

wc 

Chart  (Appendix  Fig.  A-l),  U(T.)  ,  C(T.  )  soundings  were 
adopted  as  isothermal  from  T  =  T,  (k  =  1 )  to  T  =  TQ  (k  =  0) , 
and  then  extended  from  the  last  "sounding"  point  (U[(0,  10),  T^]) 
along  the  last  emissivity  value  for  T  <  T-,  =  210  deg  k  where 
the  emissivity  of  the  water  vapor  is  known  to  be  temperature- 
dependent.   This  value  of  emissivity  is  denoted  e„  and  is 

wc 

written,  after  Sasamori  [1968] 

e^rutClO),^]  =  8.3^  T^'353  lo£l0  U(°>10>  "  W 

*  U(0,10)[-*031155  log10  U(0.10)  -  .705] 

*  [8.0/C.353  log10  U(0,10)  +  3.56)] 

+  .07262{(1.  -  .62556[U(0,10)  +  .0286]*26) 

*  [log10  C(0,10)  +  1.064]}      (A-3) 
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Figure  A-l.   Depiction  of  e    given  by  Eq.  (A-3)  on  the 
Yamato  Radiation  Chart.   Area  above  heavy  dark  curve  repre- 
sents F   +  above  level  k=l  from  a  radiative  sounding, 
wc 
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Here  the  final  term  in  the  braces  of  Eq .  A- 3  is  almost  the 
same  form  of  the  C02  correction  to  the  total  emissivity  is 
in  Eq.  A-2 .   In  Eq .  A-3,  however,  this  term  refers  to  the 
full  layer  (0,  10),  the  first  term  is  temperature  dependent 
as  well  as  dependent  upon  the  summed  absorber  masses  relative 
to  the  reference  level  k  =  10.   It  can  be  seen  clearly  that 
Eq.  (A-3)  is  usable  with  proper  final  parameters  U(0,  k) , 
C(0,  k) ,  and  T^  for  the  determination  of  the  curvilinear 
portion  of  the  downward  flux  contribution  through  level  k 
for  all  of  the  reference  level  cases  to  be  investigated  below, 
that  is  k  =  2,  6. 

It  should  also  be  noted  that  e   has  been  integrated  over 
the  temperature   T  <_   T-,  for  which  the  water  vapor  emissivity 
curves  are  temperature  dependent  (see  Appendix  Fig.  A-l)  so 
that  the  result  of  Eq .  A-3  is  really  the  average  of  the 
final  curvilinear  emissivity  value  multiplied  by  B-.  =  St  T,  . 
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APPENDIX  B:  FORMULAS  USED  IN  DEPICTING  THE  DISPOSITION  OF 
INCOMING  SOLAR  (F(A))  FROM  LEVEL  k  =  2  TO  THE 
EARTH'S  SURFACE 


These  formulas  are  depicted  by  Figs.  7,  8,  and  9  in 
Section  IV.   They  are  used  in  this  study  to  calculate  amounts 
of  insolation  absorbed  by  the  atmospheric  layers  and  the 
earth's  surface.   For  further  explanation  and  key  to  these 
formulas  consult  Section  IV. C.   Dashed  separation  lines  are 
used  to  divide  those  equations  involving  absorptions  in 
successive  layers  e.g.,  (2,*!),  (4,6),  (6,8),  etc. 


1.   Case  (1,1)   Overcast  in  Both  High  and  Low  Cloud  Layers 
(see  also  Fig.  7) 


FH  +  =  F(A)(1-.271[U(2,4)  Sec  Z]'3°3 


F4+  =  F4+(RA(1)) 

A2*H  =  F(A).271[U(2,^)  Sec  Z]'3°3 
A2*H  =  F4+.271[U(2,1D  Sec  Z]*3°3 


(B-la) 


F6   =  FlH(l-RA(l)-A(l)) 

A46  =  FlR(A(D)  (B-lb) 

A26  =  A2  4+  +  A24  +  +  A^l  6 


TD68  =  l-.271[U(6,8)5/3l'3°3 
TD910  =  l-.271[U(9,10)5/3]*3°3 


(B-2) 


where  TD  is  a  term  to  denote  the  Manabe-Moller  transmissivities 
for  insolation  beneath  a  cloud  layer.   In  this  event,  Sec  Z 


lHO 


Is  changed  to  the  mean  slant-path  sec  z  -  5/3  [after 
Katayama,  1966]. 

P84-  =  F64-(TD68) 

F8  +  =  F8l(RA(2)) 

F6+  =  FB+(TD68) 

F6++  =  F6t(RA(l)) 

F8  +  4-  =  F6l  +  (TD68)  (B-3a) 

A684-  =  F64-  -  F84- 

A68+  =  F8t  -  F6  + 

A684-+  =  F64-+  -  F8i4- 

A68  =  A684-  +  A68t  +  A684-  + 


F94  =  F8+(1-RA(2)  -  A(2))  +  F8++(l-A(2) ) 

(B-3b) 

A89  =  F8+(A(2))  +  F8++(A(2)) 


F10I  =  F9+(TD910) 
F10+  =  F10i(a(G)) 
F9+  =  F10+(TD910) 
F94-4-  =  F9+(RA(2)) 
F10++  =  F94-+(TD910)  (B-3c) 

1 

A9104-  =  F94-  -  F104- 

A910t  =  FlOt  -  F9+ 

A9104-4-  =  F9i4-  -  F104-4- 

A910  =  A9104-  +  A910+  +  A910  +  + 

A810  =  A89  +  A910 
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2.   Case  (1,0)  Upper  Overcast  Only  (see  also  Fig.  8) 


A2*H  =  P(A).271[U(2,4)  Sec  Z]*3°3 

A2Hi    =  F(A){1-.271[U(2,4)  Sec  Z]'3°3}RA(1) 

xt271[U(2,4)  Sec  Z]'3°3} 
A2*J  =  A24+  +  A24  + 


(B-Ua) 


F64-  =  F(A){1-.271[U(2,4)  Sec  Z]  *  3°3}  [1-RA  ( l)-A(l)  ] 

A46  =  F(A){1-.271[U(2,4)  Sec  Z]-3°3}A(1)  (B-Hb) 

A26  =  A24  +  AH6 


F1(H  =  F64-(TD6lO) 

F10  +  =  F10+(a(G)) 

F6  +  =  F10+(TD6lO) 

F6  +  4-    =    F6t(RA(l)) 

F1CH4-    =    F6  +  4-(TD6lO) 

A6l(H    =    F6  +    -    F10I 

A610+  =  FlOt  -  F6t 

A6104-4-  =  F6ii  -  F104-I 

A610  =  A610+  +  A610I  +  A6104-  + 


CB-Hc) 


1*12 


3.   Case  (0,1)  Low  Overcast  Only  (see  also  Fig.  9) 

A26+    =  F(A).271[U(2,6)  Sec  Z],3°3 

A684-  =  F(A).271([U(2,8)  Sec  Z]  *  3°3-[U(2 ,6)  Sec  Z]'3°3} 

F84-  =  F(A){1-.271[U(2,8)  Sec  Z]'3°3} 

F8i  =  F8+[RA(2)] 

F6+  =  F8+{1-.271[U(6,8)  Sec  Z]*3°3}  (B-5a) 

F2+  =  F8f{l-.271[U(2,8)  Sec  Z]*3°3} 

A68+  =  F8+  -  F6+ 

A26l  =  F6+  -  F2+ 

A26  =  A26I  +  A26+ 


F9   =  F84-[1-RA(2)  -  A(2)] 
A89  =  F8+[A(2)] 


(B-5b) 


F104-  =  F9+(TD910) 

FlOt  =  F10+(a(G)) 

F9+  =  F10+(TD910) 

F9+4-  =  F9  +  (RA(2)) 

F104-4-  =  F9  +  MTD910) 

A9104-  =  F9+  -  F104- 

A910+  =  FlOi  -  F9+ 

A9104-4-  =  F9  +  I  -  F10+  + 

A910  =  A9104-  +  A910i  +  A9104-  + 

A610  =  A68I  +  A68+  +  A89  +  A910 


(B-5c) 
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